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OPTIMIZ/TIOK HT Mlffil-SPIITOIE /IJTOM/iTICS 
- A PB0B/3ILISTIC APPBDACH 

An analytical methodology is proposed to determine 
optimum cutting conditions and optimum tool-replacement proce- 
dure for manufacturing operation performed on multi-spindle 
automatics considering tool-life and cutting-force as proba- 
bilistic in nature. The optimization is carried out w.r. t* 
either minimization of cycle time or minimization of total cost 
of production^ considering the various constraints imposed by 
machine and cutting-tool characteristics, physical process, and 
the job specifications. The concept of diance constrained 
programming has been utilized to account for the probabilistic 
nature of tool-life and cutting-force, 
r 

\ The model ooQisiste of a non-linear objective function 
subject to several non-linear constraints. Ai interior penalty 
function ^proach is employed to help use Sequential Unconstrained 
Minimization Technique (SDMT), Davidon - Fletcher - Powell 
Variable Metric Method (itFP) is used to advaitage with Golden 
Search Technique for carrying out unidirectional minimi zat ion Tj 
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The piDpcsed raeiiiodology is computerized to deter- 
mine optimum, machining parameters aid to detemine the total 
cost of production for the tv/o tool- replacement procedures. 

The piocedures are : l) unscheduled tool- replacement aid 2) 
scheduled tool- replacement , It is shown that the optimum 
cutting conditions, cycle time, and the total cost of pixiduc- 
tion are significantly effected by the probabilistic nature 
of coefficients in the constraints. The computer package 
has been tested through an example of a job to be produced 
on multi-spindle automatic. 

In this study the effect of l) capacity of machine 
tool and 2 ) the type of machine tool depending on T/fcether 
all the spindles have common speed and feed or the spindles 
having different speeds aid feeds has also been investigated. 
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UTROnJCTIOI 

The concept of manufacturing is undergoing a revolution. 
The- traditional thinking that maauf actu ring is aa art is, slowly 
but surely, being replaced by the idea that it is a scientific 
branch of engineering. To quote Herdiant (33)5 

"Basic and long-rmge technical changes are today 
challenging the conventional economics of manufacturing. 
The key weapon to vhida industry is turning to combat 
this challenge is new optimization technologyi computejc- 
related technology and versatile automation are its 
main components. The prime technical advance vh-ich has 
made possible extensive use of both of these categories 
of nev; technology is the digital computer; manufacturing 
is apparently already its second largest area of appli- 
cation in industry.’ Versatile (variable programme) 
automation, under the impetus of computer’s capabilities, 
is now striving to upgrade its cone rol capabilities from 
tho se afforded by feedback and supervisory control to 
those offered by self - optimizing control. Simulta- 
neously, it is rapidly upgrading its programme input 
sources from previously restrictive manual and analogue 
inputs to digital inputs. This process is leading to 
the ultimate tool for optimization o f manufacturing, 
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direct digital control, in vdiida the control loop is 
closed through the computerj the fully self - optimizing, 
fully aitQPigted miSaufacthring system is its pot^tial 
climax . . , 

Discrete item mahttfacturing li*s laag been re’CeptisEc to- 
automation innovations. Since 1955, one of the fastest growing 
areas is the . development of Computer Assisted Manufacturing (OAM) 
systems. The control of digital computer has been fully accepted 
in the continuous or process industries. However, this acceptance 
has not heen so prevelant in the discrete item industries for cer- 
tain rather obvious and complex reasons. 

An optimizing machine can be thought as a separate sec- 
tion in a fully automated system. The purpose of the machine is to 
optimize the manufacturing activity in its entirety. The optimizing 
section may further be divided in order to handle various aspects of 
manufacturing - raw material control, production planning, schedu- 
ling, production, handling of finished products and so forth. The 
optimizing equipment must be equipped with certain logic, ana- 
lytical or heuristic, in oi^er to arrive at optimum decisions. 

The present vroik deals with the manufacturing aspect of 
a concern producing discrete components. The ability of any manu- 
facturing firm to successfully compete in the market is highly 
dependent on the cost of its manufacturing processes. To quote 
Cook ( 9 ), producing satisfactory parts at the lowest 



possible cost cai be called the 'First law of Production' " 

In most oiganizations those decisions which directly relate to the 
cost of producing a product are made on the basis of experience 
(often limited) and ?ri.th the assistance of various charts and tables 
of "proper operating conditions" put out by the manufacturers of 
machine tools, cutting -tools, etc. Although the basis for rational 
selection of manufacturing parameters was set down by I ay lor in 
the early 1900 ' s, only recently -people have attempted serious 
utilization of analytical selection methods. 

Considerable woik has been done applying analytical and 
heuristic approaches for single - tool operation using determinis- 
tool-life values. It has been recommended that tool - life predic- 
tions be made on a probabilistic basis. However, the economics of 
multi - tool operation -with probabilistic tool - life consideration 
has received vezj'' little attention, possibly because of its inherent 
complexities. It was thus felt that stiucturing of a model for 
multi - tool operation and subsequent ^plication of a mathematical 
progranming technique would be a worth-viiile effort. Therefore, 
the present vroik is directed towards the optimal selection of cutting' 
parameters for multi - spindle automatics when the tool - life is 
considered probabilistic. Further, the woric includes the selection 
of tool«#replacement programne for changing the cutting -tools. 

The problem of optimization for multi - spindle automatics 
is extremely complex due to diversities in the type of operations 
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performed at each station. This problem is further aggravated due 
to the probabilistic tool failure characteristics of the tools 
involved. The problem as such, consists of having a methodology 
to provide answer for the following decisions ; 

1. Determination of machining parameters for 
each station, 

2. Determination of optimal tool - replacement 
policy, and 

3. Selection of a machine tool. 

To seek a logical answer to the above problem, an analy- 
tical model has been formulated. The stochastic nature of constants 
and exponents in tool - life and cutting - force equation is accoun- 
ted for by structuring the problem as a chance constrained progra- 
mming problem. The constraints are imposed due to machine and 
cutting - tool characteristics and job requirements. The objective 
functions studied are : 

1. Minimization of cycle time, and 

2, Minimization of total production cost for a given 
production lot size. 

Ohe model consists of a non - linear objective function 
subject to non - linear constraints. M interior penalty hinction 
approach is employed to help use Sequential Unconstrained limimiza- 
tion Technique (SUIT), Davidon - Pletcher - Powell Variable Metric 
(dT!P) method is used to advant^e wilh Golden Search Technique for 
carrying out unidirectional minimization. 
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The algorithm has been piogrammed in. iOHTEiW IV language 
for use in IBM 7044 digital computers. The computer package is 
quite general in nature to handle optimization problem for manufac- 
turing operation in vihich material is removed in form of chips. 
However, the package is tested only for a job to be produced on 
multi - spindle automatic machine tool. 

The present work is but a step or two leading to the final 
goal viz. fully automated manufacturing i^stem. However, 1he even- 
tual realization of such a system will and must be an evolutionary 
process. 
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EITEEAEUEB SUEW 

The art of metal remo-ving is as old as human ci-vilization 
itself, Kevi? chapters are being added continuously to meet the 
challenges of time. The state of the art that is mtnessed today 
is constituted of a number of broad areas each of them forming a 
separate field of its om. The areas relevant to the present work 
are s 

1, Tool - life and its Characteristics, and 

2, Economics of Machining, 

The following sections discuss the literature regarding 
the above mentioned areas, 

2,1 Tool - life and its Characteristics 

T'. ¥, Taylor (43) was the first to investigate the tool - 
life relationships. The well known equation relating tool - life 
and cutting - speed does not, however, take into accoimt the magni- 
tude of the chip cross - sectional area and as such a different 
’Taylor Constant' is required for each chip cross - sectional area. 
The necessary refinement in Taylor's equation has been done by 
Kronenberg (30) , He derived initially an elementary and subsequently 
an extended cutting - speed law. The extended cutting - speed law 
considers the effect of both chip cross,- sectionaJ. area aid 
slenderness ratio, -vifoile the elementary law does not include the 
effect of the latter. Kronenberg, however, pointed out that chip 
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cioss-sectional airea has a considerable laiger effect on the tool 
life than the slendemess ratio. 

In a real-life situations tool-life rarely coincides 
with the predicted value. Yet, despite practical data on toolr 
life values, veiy little has been known about actual magnitude of 
•tile scatter. Various explainations have been forwarded to explain 
this scatter. Eimer and WU (l5, 16) have attributed "the variance 
in tool-life purely to the expeilmaatal error. On "the o'tiier hand, 
s'tudies conducted by Wager and Barash (46) have led to the conclu- 
sion -that -fche variation in tool-life can not be attributed to 
'esperimaital error' but it is "the resultent effect ofmgjay factors, 
material hardness, workpiece rigidity of mounting, themo-electric 
effect, etc,, Al-fchough -the effect of each individual factor is 
negligible, the sum total assumes appreciable proportion. The 
B'tudy also reveals that -(he variations in tool-life are inherent 
to the physical nature of the process vhidi, like so mmy oiher 
physical processes, is stochastic. The aithors analysed experi- 
mental data and found 'that the life of High Speed Steel (HSS) "tools 
would be best represented by a positively skewed normal distri1xition% 

Exrtensive efforts ( 17) have been directed to determine 
tool-life data for single-point tool but the available literature 
reveals that the tools of complex geometry (e.g., form tools) 
have recei'ved very little attention. 
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2.2 Economics of Machining 

The Vv'ork of Taylor (43) in the United States, culminating 
in a major paper in 1906, had aa important Influence on the develop- 
ment of economics of machining. He suggested that a relationship 
must exist between machining time and tool - life in order to have 
minimum cost of production. Till to - date this particular facet 
of metal removing has received concerted effort to establish con- 
ditions with a view to manufacture economically. 

The discussion on the economics of machining is presented 
in a hierarchical order. 

2,2.1 Single - Point Tool Economics 

Many authors have perfoimed experimental and analytical 
studies of the single point tool case to understand better how 
the cutting conditions affect production times and tool lives, and 
hence the costs of cutting. The problem consists in detennining 
machining parameters in order to optimize a given objective funo- 
tion. The variety of objective functions in common vogue are : 

1. minimm cost per piece, 

2. minimum machining time per piece, 

3. maximum profit per piece, 

4. maximum profit rate, and 

5. meeting of demand aiccessfUlly, etc. . 

The flood of Uteratore in this field can be classified 
under two major sub-headings depending upon the type of formulation 
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and methodology used : 

(a) Unconstrained Optimization, ;>nd 

(b) Constrained Optimization. 

(a) Unconstrained Optimization J 

This optimization methodology has been adopted by a number 
of researchers including Taylor (45a), Brewer (7), Gilbert (22), 
Weill ( 47 ), Gilman (23) and Golding (s). They utilized empirical 
fomulae to express tool - life as a function of various machining 
parameters such as cutting - speed, feed aid depth of cut. The 
objective function is differentiated to obtain the required optima. 
These models aire quite introductoiy in nature. However, these ini- 
tial efforts triggered a lot of interest among the research workers 
to develop improved models. 

Held, et al (l9, 20) have taken tool - life data add 
have used a digital computer to generate tool - life equations. 
These equations were used to develop relationships for the determi- 
. nation of optimal speed. Only the first three criteria have been 
used in their analysis . 

Wu & Exmer (48) concemed themselves with deteimining 
speed and feed so as to maximize profit. This was accomplished 
by assuming functions for marginal revenue and marginal cost, 

!3^he optimum machining paraaeters were calculated by equating the 
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expressions. Armarego and Eussel (l, 2) have argued that for a 
given income per component, profit ma(?e by the manufacturer per 
unit time is not optimum #ien the criterion of maximization of 
profit is used. Hence, from manufacturer’s point of view without 
increasing the cost of commodity is of prime importance. A similar 
analysis has beem made by Okushima and Hitomi (35). Wu and Tee (SO) 
have investigated the profit concept to obtain a better understand- 
ing of the profit region. The authors have found that optimum 
cutting speed seems to lie between minimum cost speed and the speed 
for meeting' the demand satisfactorily, 

Emer and Eaiia - Gonzalez (l4) have studied the sensiti- 
vity of machining cost with respect to the cutting conditions. This 
study was primarily conducted to understand the shape of the cost 
curve in the neigh bouiiiood of the optimum, cutting conditions, 

in order to account for tool life scatter, Ermer and lu 
(l6) suggested that optimal cutting speed should not be defined 
uniquely aid instead a range daould be associated with it. They 
advocate the minimax principle for the purpose of speed aid feed. 
Elmer aad Morris (l5) in another paper have recommended the use 
of a correction factor to account for uncertainties in tool - life. 
They claimed that the application of correction factor is more 
direct and simpler con^jaxed to the minimax piinciple, 

Vul’fson and Deiyabin (44) developed an iterative pro- 
cedure for selecting the speed and feed viiich yields minimum 
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cuttizLg cost. A computer based procedure was designed specifically 
for numerically controlled lathes, 

Yeiy few research workers have directed their efforts 
to machining operations o1her than simple tuining, Aimarego and 
Hissel (2) has presented an analysis for single - pass taping and 
milling processes based on maximum profit criterion. 

b) Constrained Optimization • 

In a real life situation the choice of machining para- 
meters is constrained by machine and cutting - tool characteristics 
and "the component specifications. Hence, a constrained problem 
would represent machining activity in a more realistic way. 

A linear programming formulation has been developed for 
the purpose of selection of a machine tool and corresponding cutting 
conditions such that a rotational part may be manufactured with the 
least cost in a job - shop environment (js). Krishna and Berra (28) 
has formulated a multi - pass face milling optimization model so 
as to minimize the cost of production. The model consists of a 
non - linear convex objective function constrained by linearized 
constraints. Bosen’ s G-radient Pro jection method is used to solve 
the problem. 

Hiattachaiya and Paria - Gonzalez (5) have presented a 
model considering airface fini^ as a principal constraining para- 
meter. In their analysis the slack variables are introduced in 
the power form and iiie method of Lggrange Multiplier is used for 
ascertaining local minima. 



In 1972 Iwata (25) proposed an analytical method applying 
a chance constrained programming concept to determine the optimum 
cmditions introducing stochastic concept in both the oboective 
function and the constraints. He has advocated that the constants 
and exponents in tool - life and cutting - force equation daould be 
considered as stochastic paraneters. It is ^owi that the optimum 
cutting conditions are significantly affected by these considera~ 
tions . 

2,2»2 Economics of Multi - Tool Set ups 

This area remains, relatively, unexplored because of a 
few additional problems associated with it. One of. the major 
problems in multi - tool set ~ ups is the identification of opti- 
mum groups of operations. The problem of grouping is irrelevant in 
the case of single-'tool . The over all problem of nulti - tool set 
up can also be divided into the folloiying sub - problems ; 

a. Identification of optimum groups, 

b. Determination of optimum cutting parameters, and 

c. Detennination of tool - replacement schedule. 

The sub' - problems are discussed in detail in the follow- 
ing sections. 

a) Identification of Optimum Groups 

A paper on the related subject of determining groups 
among the operations to be performed on a particular job is due 
to Jones and Morgan (27). They have employed a grajhical technique 
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called Analysis of Interconnected Decision Areas originated by 
luclanan and Stringer (26). The method has no objective function 
and thus does not seek optimality , but rather it attempts to 
identify feasible alternatives when they exist. 

Bartalucci (3) has formulated a dynamic programming 
model to solve the problemj but as number of operations increase 
the numerical effort required by the method increases enormously. 

b) Determination of Optimum Cutting Paremeters 

Brown (7a) revieyred the literature and opined that -vdaen 
an operation comprises of two or more passes, all passes should be 
of Same depth so as to achieve an optimum solution. 

McCullou^ (32) has reported- a very elouentary and 
approximate analysis for nulti - tool case. The equations derived 
were said to express the optimum tool - life for minimum cost and 
maximum production rate. 

A detailed procedure to find optimum cutting conditions 
for multi - tool operation is presented by Goranskii (24) in his 
book ’Theory of Automabion of Production Planning and of Tooling' . 

A major contribution of Goranskii lies in linearizing the model 
by taking logarithm of the both objective function and the cons- 
traints. To summarize, the algorilhms put forward suffer from ■ 
two major drawbacks. First, the tool - life equation has four 
constants and eight exponents vhidh in real life situ^ion are veiy 
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difficult to determine. The second reason for ccmplexily of the 
algorithm is that the linear - progi-am as initially developed, 
be infeasible and the algorithms are required to adjust systemati- 
cally to the various parameters until a feasible solution is 
obtained. 

c) Determination of Tool Replacement Schedule 

In machining operations involving many tools simultaneous- 
ly, such' as in multi - spindle machine tools, the choice of tool- 
change procedure and tool*change time interval is critical to the 
overall efficienqjr and economy of production. The study of the 
literature available reveals that the problem of determination of 
optimal cycle - length of a cutting-tool presumes machining para- 
meters to be lonovsn. 

A cost model for strict interval - replacement policy 
has been suggested by T aha (42). The objective function includes 
the cost of defective parts produced in addition to the tool - 
replacement costs. In order to simplify the model it has been 
assumed that any increase in percentage of defective items, over 
the situation when 'the process is stable, is caused only by a poor 
condition of the tools. 

Ebert and Hershauer ( 12 ) has conducted a simulation 
experiment to examine cutting - tool operating costs under vary- 
ing conditions. In their model the tool-replacement policies 
incorporate." different type of (jxality control decision rules. 
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The strict interval ~ replacement policy has been found signifi- 
cantly better than the policies governed by quality ccntrol 
decision rules, 

Okushima and iUjii (34) has formulated cost models for 

* 

three different types of tool - replacement schedules. The author 
has recommended that due to inherent uncertainties in the tool- 
life, the optimum solution can not be reached directly. Hence, a 
simulation method offers a useful means to determine an economical 
tool - daange sciiedule. The computed data indicated that the 
scheduled tool - changes are more economical than the policy of 
replacing the tools as 1hey fail. 

luncan (ll) has developed a very sophisticated mathema- 
tical and numerical model for scheduling cutting tool - changes, 

2,2,3 Grroup Technology and its Bamifications 

G-roup technology is a technique of grouping similar 
components together in order to eliminate repeated efforts re- 
quired to design and plan for a new component vhich resembles 
another component already designed. 

Opitz ( 36 ) has developed a coding procedure for the 
jobs based upon the philosophy of part family manufacturing, 
production method that involves machining of parts in fanilies. 


* Explained in section 3.1»6 of Chapter III 
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Scott ( 41 ) has reported, a method of automated pleaaning for manu- 
facturing which is called 'Regenerative Shop Planning'. This 
system put an end to the repetitive manufacturing decisions, 
redundant operation plannitog aid the piling up of paper woifc. 

In this system the logic by wiiich the original component was 
manufactured is placed in computer storage and utilized in plann- 
ing of subsequent components. In order to institute this method 
■the aithor has stated that one must be able to generate part 
families which have many characteristics in common. 

As recently as August '73 Eajagopalan (39) has intro- 
duced a graph theory approach to the problem of grouping the 
machines in a batch - production shop into cells in such a way 
that each cell takes care of a family of similar components. 

2.2.4 Automated Manufacturing Systems 

Most of the work done in this area is of proprietary 
inha'ture, therefore, has not been reported in literataie. Some 
of the well publicieed systems ©re APT ( 4 ), MII&iAP (4), AUTOPIT ( 4 ), 
EXAPT II ( 4)5 etc.. These systoas allow the description of a 
part to be made in a language that closely resembles english. 

Berra and Barash, and Krishna (4, 29) have brought forward a com- 
prehensive system for the automated plaining and optimization of 
a manufacturing process. Their interest has centered on numeri- 
cally controlled lathes, but the scope of their system ranges 
from 'the specification of the component to 'tiie numerical control 
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of manufacturing operation. Their -work has laid foundation for 
automated planning but laclcsin generality. The essential reason 
is that a numerically controlled lathe usually has a single tool 
^lich makes a series of cuts to generate a part, -viiile multi ~ 
tool operation has got different characteristics of its own. 
Obviously single - tool aspects are embedded in multi - tool ope- 
raticn , 

To summarize, the literature survey reveals the lack 
of scientific literature in the area of multi - tool operation. 

It also opens up a new awareness of the distance between the 
ultimate and the present. 



GHJiMEE III 


SENERAL DESCBIPTIOF OE MULTI-TOOE 
OPEELAIIOE SISSmSS 

Before delving into the problem foimulation pioper, it 
is essential to have a better understanding of 1he i^stam under 
study So as to help modelling. Hence, a description of some of 
the important elements synthesizing the ^stem is delineated, 

^Qle relevant portions of iiie model involving deteimination of 
cutting conditions, cutting time and panduction cost is presen** 
ted in CJh^ter 17. 

3« 1 Elanents 

The major ^stem elanents of the production activity 

are : 

1, The machine tool - a multi-spindle automatic, 

2, The cutting tools, 

5, The woiispieces, 

4. The machine operator, 

5. The tool-setters and tool-setting activity, ' 

6. The tool re^arpening acti-';i1y, 

7. The tool failure, and 

8. The tool changing activity* 

Most of Ihem are in common vogue aid need ro esplar 
nations. Discussion only whenever pertinent has bean included. 
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5»1.1 The Machine Tool — a Multi— Spindle iiutomatic 

The machining methods employed on multi - spindle automa- 
tics depend upon liieir principle of operation and on the type of 
blank used. Depending on the type of stock or blank, the machines 
can, be classified as 'bar type' or 'chucking machines'. According 
to the principle of operation, the automatics are classified as 
parallel and picgressive action machines. In the parallel action 
sutomatics tiae sane operation is performed simultaneously on all the 
spindles. Thus, during one operating cycle as many workpieces are 
produced as there are spindles. In the progressive action automatic, 
each spindle with its blank or stock is indexed to the next position 
after machining is fini^ed at each position. Every workpiece has 
to pass throi;igh each position consecutively. 

In multi-spindle automatic the cycle time is equal to -Hae 
maximum of the minimum possible production time for each spindle. 

In such a situation the most expedient producing method is one in 
vhich the time required at each station is s®e. 

Most models of multi - spindle automatics have common speed 

and feed for all the spindles. But the modem trend is to have 

different speeds and feeds for each spindle. In the present work 

this type of automatics are categorized as type - A, The previous 

type of automatics are categorized as type - B,- 

!= machining time for the station requiring ma 
cutting time + tool approach time + indexin 
+ tool withdrawal time + loading and imload 
time for the job. 


*Cycle time 
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3*1.2 The Machine Operator 

The span of activity of the operator, generally, varies 
with aiop conditions. It is not unusual to assign more than one 
machine to one operator in case of bar type automatics. The number 
of aaciiines assigned to an operator is a function of : variety of 
activities that he ha^ to perform, the expected frequency of machi- 
ne requiring attention either due to tool failure or exhaustion of 
ravY material, etc.. 

3*1.3 The Tool setters and Tool Setting Activity 

This activity is critical in multi - tool operation. The 
practice of preset tooling is usually adopted to reduce down time. 
According to the piedetermined replacement schedule the tool setter 
replenishes the-, machine with tools. During the setting time the 
production has to be stopped, so the downtime is strong function of 
the policy regarding the replacement schedule. 

3,1.4 The Tool Re^arpening Activity 

The tools which are regrindable are ^arpened to mdee them 
u suable once again* The inventory of cutting -tool includes three 
echelons j the tool inventoried in the expense crib, (the new tool 
inventory and resharpened inventory), the tools inventoried at the 
machines, (machine float), and the tools stored in the grinding 
^op awaiting regrond, (regrind inventory). 
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Once the grind inspector receives a. shipment of worn 
and danaged tools, he has to make several decisions. First, he 
decides vtiich tools are to be reground and which are to be scrapped. 
Second, he m^ detenaine how much stock is to be ground from each 
tool. Finally, he dete mines whether or not he quality of regtound 
tools is satisfactoiy. The grinders usually set their machines up 
so that one set-up will suffice for an entire lot. This means that 
the amount of stock removal for each lot depends upon Ihe worst tool 
and so is the time of grinding. Many tools such as taps, small 
drills, throw away - tips, etc, are never reground. 

3.1.5 The Tool Failure 

The tool failure phenomena has been recognized as a pro- 
babilistic event. Inspite of the uncertainties involved, empirical 
relations have been proposed to find out expected -values of tool 
life. 

In this work, the production life of a cutting - tool is 
defined as the number of satisfactory workpieces produced from ihe 
time of tool - replacement to failure. It needs to be pointed out 
that whenever the tool is replaced before its failure, it is because 
of particular tool-replaCQuent schedule. It is called premature 
removal. It will be discussed later in the section 3.1.6 that in. 
some of the tool - replacement schedules premature removal is 


inevit able 
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Duncan (ll) has reported that the tool failure can be 
attributed to the agi^ig effects; cominonly categorised as weaTj dete- 
rioration, fatigue, etc,. 

lot of work is reported on the study of tool failure pheno- 
mena. Duncan has listed reasons for failure under two separate 
headings, 

i) Sources dependent on tool, and 

ii) Sources independent of tool. 

Delving into all the pros aid cons of the phenomena is not 
within the scope of the present research. It suffices to state that 
all the reasons grouped together make tool failure a random phenomena* 

3 . 1.6 The Tool Changing Activity 

Tool changing is primaJd-ly a corrective action to maintain 
quality and avoid accidents due to loading of wom out tools. To 
take 1iiis action, the management must have proper guidelines for an 
economic tool - replacement schedule. As per Gfcudiima and Dujii (34) 
tool-daange procedures can be divided into the following categories^ 

1 . Unscheduled tool — change : In this case, any indivi- 
dual tool vd.ll be changed whenever it fails. Sche- 
matically shown in Pig, 1, 

2 . Scheduled tool - change ; In this case, there is a 
predetermined tool - change interval, but some un- 
expected tool failui>es may require extra tool changes. 
This alternative consists of two sub-alt ematives. 
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First, all tools ’.'ri.ll be changed without regard to the 
histoiy of eadi tool during the pieced drg interval. This procedure 
is shown schematically in Fig. 2. 

Second, only selected tools will be cfianged, that is, 
those tools which have been changed during tiie preceding interval 
will not be changed at this time. Only those tools vhidi have not 
been changed will be replaced. This procedure is shovn schemati- 
cally in Fig. 3. 

In the present vnrk the two sub- alternatives are named 
as schedule 2-A and schedule 2-B respectively. In Fig. 1, Fig. 2, 
and Pig. 3, the example of machine tool witli three spindles is 
given. The white circle in the figuios represent the changes 
occuring due to tool failure and the crosses represent the tool- 
change at the predetermined tool-change time. Here, it can be 
observed that in case of scheduled tool-change procedure, the 
premature tool removal may occur quite often. 


I W , "V- . . ^ P : — ^ — ’■ — ■ 

Unscheduled lool-Eeplacement Scheduled Tool-Replacement 2-A 

FIG. 1 FIG. 2 



Scheduled Tool -Replacement 2-B 
FIG. 3 
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iny tool clisngiGg activity is associated with downtime 
production losses. Often, the costs associated with these down- 
times are intaigible in nature. The time rtcuired for tool chaiging 
noimally depends only upon the relative position of the tool on 
slides, tool inventories aad quality control tests after the setting 
of the tool. The change - over time can be split into j 

Time for, 

(1) machine shut down, 

(2) tool retrieval, 

(3) tool removal and replacement, and 

(4) quality control testing. 



CHAPTER 11 


REDBIEM EOEMnLATIOI 

4.1 The Pixjblem 

It is req.uired to have a melhod of deteiminiag an optimal 
strategy for the machining operations performed on multi-spindle 
automatics. This involves the machine tool selection, tool-replace- 
ment schedule and the determination of machining paraaeters - speed 
and feed rate for every tool. The objective is to minimize either 
the total cost of production for a given lot size or total machining 
time • 

The problem requires the following types of data as input 

1, Job specification as specified by the customer, 

2. Operation plan of the job to be manufactured. The operation 
plan specify - (a) groups among the operations, and (b) assign- 
ment of spindles to the various groups formed. 

5. The relevant data generated from job drawing and its specifi- 
cations, 

4. The data a^ou-'t cutting - tools used include, 

(a) the type of the tool used, 

(b) the mean and standard deviation of the tool - life 
and force equations exponents for each tool under 
consideration, and 

(c) the relevant characteristics of the tools used. 
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5« Machine type and its characteristicSj 
6, Cost and time parameters, and 

I. Size of the lot to be produced and the attrition allowance, 

Clicptcr VI contains a detailed listing of various 
characteristics for all types of data mentioned abo-ve. 

The following sections provide a comprehensive descrip- 
tion of the various irodels formulated to solve the problem. The 
formulation is based on the assumption that the optimization is 
carried out only for primary operations like, turning, drilling, 
form tuming, parting off, etc. The operations like, thread cutting, 
chanfering, knurling, reaming, etc. are taken as secondary operation. 
These operations, if at all to be performed on a job, are excluded 
from the optimization routine. 

Separate models are formulated and studied two types 
of aitomatics - (i) multi - spindle artomatic for viiich different 
speeds and feeds are available for each spindle termed as Type-A, md 
(ii) the machine tool having common speed and feed value *r all 
the spindles termed as Type-B, 

4,2 The Objective i\mction 

The objective functions for optimization of machining 
operation in common vogue are : 

1, Minimization of cost per unit, 

2^ Minimization of production time perimit, 

3» Maximization of profit per xmit, 
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4. Maximization of profit rate, and 

5. Meeting the daaand successfully. 

The profit finction can not be justified becaise of its 
subjective nature, A distinct estimation of contribution of each 
machining operation to the total profit earned is very difficult. 
Hence, the decisions based on the profit function may invol've hi^ 
degree of arbitraryness and the viiole optimization efforts loose 
their value. 

The demand meeting function may be appropriate only for 
a particular shop - condition : a shop may have to pay a very high 
penalty for making delays in the ^ipment. 

Therefore, ihe present v/oik restricts itself to study the 
follovving two objective functions. The objective functions are ; 

1, Minimization of total maciiining time. In the multi - spindle 
automatic operation the cycle time dictates the total machining 
time. Therefore, minimization of machining time is equivalent 
to minimization of cycle time, aid 

2, Minimization of cost of production for a given lot size. 

4.2.1 Cycle Time Minimization 

The various components constituting cycle time can be 
listed as follows ; 

1. maximum of machining times for all the cutting 
tools, minutes, 
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2. t - time required for withdrawal of all the tools, 

minutes^ 

3. tj - time required in. indexing of the spindles, minutes^ 

4. tj^ - time required for unloading and loading, minutes, 

and 

5. t - time required for approach of the tools, minutes. 

The cycle time, TC, in minutes can be expressed as 

TG = t. + t + max (t. .) + t + t (4-.l} 

1 a . . 13 w I 

2-?3 

i="]52 ,.«. 5 i^ 

3 == 1, 2,... 5 n^ 


#iere 


i - index for spindle, 


3 - operation index for operation on i spindle, 
t. . -• machining time of an operation, minutes, 

13 

II - number of spindles, and 

• th . 

n. *• total number of operations performed on i spindle* 

As it is pointed out above that t^j t^^ t^^ aad t^ time 
elonents are constant dependiag upon a combination of machine tool 
and the job to be pro duced^ hence need not be considered in the 

objective function. Ihus the objective function reduces to, 

^ max “ - ) , t ^ i j 2 j * # * , 

i,3 

3 — 1 j 2, •« » , h^ 


minimize TC 


(4.2) 
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In liie following sections the cycle time will be understood 
equal to its variable part onlj''. 

The machining time of an operation can be expressed in 


terns of machining parameters as follows 


t. . 


V. 

1 



(4.3) 


viiere 

"bl3- 

1. . - tool travel length for j operation on i 

spindle, inches 

~ i”^^ spindle speed, rpm., snd 

^ j. .th . „ .th 

f . . , “ feed rate on i spindle for j operation, 

X J 

indh/rev. 


It is to be noticed that the subscripts for speed and feed 
rate are meaningful only in the case of automatic type - A, but the 
subscripts do not affect the fomulation for the automatic of type -- 
B. The affect of machine type will be realized only in the develop- 
ment of constraint set, 

4,2.2 Production Cost Minimization 

The major components contributing to the total cost of 
machining operation are ; machine and operating cost, tool costs 
and non productive costs. The machining cost is that of operator 
and overhead computed over the production time for a lot to be 
produced. The tool cost is comprised of tool edge cost, machine 
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setting cost 5 intemediete tool - changing cost, and tool regrinding 
cost, Pinally, the non productive cost includes, the production 
lost in non-cutting time due to intermediate settings because of tool 
failure and the cost, of operator for the same period. 

In the following section the expression for various costs 
are developed for the two tool repla,cement procedures; l) unscheduled 
tool - replacement (Pig. l). Case 1, 2) scheduled tool-replaoement (2~A) 
(Pig,2) » Case 2 • In Case - 1 there will be unexpected tool changes 
and may become expensive and time consuming in comparison to schedu- 
led ones. 'But excessively frequent scheduled changes are also un- 
economical. 

Case 1 - Unscheduled lool-Replacement : 

Machining and Operating Cost 

The total machining and operating cost, C^, in rupees can 
be expressed as s 

C = TO . . (C, + C„) (4.4) 

m I 

vdiere ' 

KJ - total number of jobs to be produced (the lot 
includes the attrition allowance)^ 

- direct labour cost, Rs./min., and 
Cg - overhead rate, Rs./min.. 

Tool Cost 

The uncertainties involved in tool - life prohibits one 
to know in advance 1he moments of tool failures. In such a case 
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expected values can be priced out only. The expected value of tool - 
life depends upon machining paraieters. Hence, tiie selection of 
machining pajrameters should prolong the tool - failure and reduce 
tool replacement costs. 


Vyoskovoskii (45) has reported that the case of unscheduled 
replacement is seme as series replacement, A system comprised of 
many components of average finite life has an average life vhose 
reciprocal is equal to the sum of reciprocals of life of individual 
components. In the case of multi - spindle automatic the machine tool 
can be thought of as a system and the components of the system are 
various cutting - toolshaving a finite mean life. In the present 
TOrkjthe machine tool failure is ignored. Hence, the average life, 

T, of the system can be expressed as 


J. 



■^1 

D=1 



(4.5) 


where 

“fch. 1 j1i 

T . . - mean life of j cutting tool on the i spindle, 

A cost expression for the system is next formulated. The 
formulation has been carried out with the assumption that ihe various 
downtimes enumerated in . the section 3. 1.6 of Chapter III are combi- 
ned into are single term - tool changing time. 

The expression for the total tool cost, CGH, is writtsn 


as 
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n. 




n. 


CCH 


i=1 D=1 




m. . + 1 
ID 


+ 1 


EG. . + 
ID 


x (r 

i=1'- D=1 ' 


p . . + 

ID 


1 

1 


.ri 

m + 1 • ! . 

- Ji 


. TG. . . CG + ? p. . I , TS. . . GS + TS . C 


ID 


L. 


ID 


ID 


(4.6) 


■wile re 


IP. 


PiD 


JJl 


I. . 
ID 


GG ~ grinding labour rate, Rs./min,, 

IS - initial machine setting time, min,, 

G^ - machine tool setter rate, Rs,/min., and 


- indicates that greatest integer value of Z is to be taken* 


"feii 

lor the i spindle i 


.th 


I. . , UJ 


IP. . “ total actual machining time for j tool - - - »- , 

XD Vi . fi . 

1.. - actual machining length vdth tool j, inches, 

X J 


mm. 


m. . 
XD 


number of possible regrinds mth tool 
-th 


EG.. - cost of anew j tool - edge, Rs. , 

IG.. -- grinding time for a cutting tool j, minutes, and 
IS. . - inteimediabe tool - change time for 3 tool, min,, 

j- J 


Eon-P reductive Cost 


The non-productive loss occurs vhenever the machine is stop- 
ped either for the purpose of tool - change or because of machine 
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maintenance. In literature it has been found that this cost is 
accounted in teims of direct labour w^e rate. However; the loss 
incurred is not only the money paid to the labour but there is a 
production loss also during the machine stoppages. The loss of 
production can be accounted by imagining of the possible output 
during the idle time. Therefore; the loss of production due to the 
down time can be considered as lost profit. Since, this factor is 
highly nebulous and the formulation assumes that this information 
will be provided by the management. Hence, the non - productive 
cost, CIP, is expressed as follows } 

TS. . 

( 

i=1 3=1 


H n. 

CMP = ( 


JaL 

TO 


^ • .QiS J 


OP 


M 


1=1 0=1 


n. 


i3 ! 


TS. . 
10 


(4.7) 


where 

CP - expected profit per unit, Rs, 

The total cost of production, CT, is summation of equations 
(4,4), (4.6), and (4,7), Therefore, 

CT = C + GCE + CMP (4,8) 

Case 2 - Scheduled Tool - Beplacanent 

Scheduled tool - replacement is characterized by pre- 
determined time interval after which all the tools are replaced 
irrespective of their condition. Moreover, if any tool fails before 
the predetermined tool - replacement interval, it is also replaced 
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at the time of its failure. The knovdedge of probability distribu- 
tion function (p.d,f.) of tool-life cai help in the detenaination 
of total number of replacements before the scheduled time of tool- 
change, But veiy little experimental vroik has been reported in 
literature to establi^ the exact nature of p.d.f. of the tool-life. 
However, the p.d.f, will be a composite distribution obtained by 
taking weighted mean of the distribution of different nodes of tool 
failure 5 wei^tages being assigned on the basis of relative dominance 
of each mode of failure. The unavailability of tool-life distribu- 
tion for various types of failures compels one to adopt a simple but 
approximate method for the purpose. In this study a factor for 
tool on i"^ spindle is defined by the following equation ! 



( 4 . 9 ) 


where 

-2 ^ 9 < 2 

SI - tool-replacement interval, number of componoits 

produced during liie interval and 

- "fell 

^T. - standard deviation of tool-life of j tool on 

i”*^ spindle . 

There are three possible values of C. .’s, -.They are ; 

1. 0 ^ C. . 4 1 , . 

2 . 0 ^^ > 1 » 


3 


C. . < 0 

X3 
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The first possibility does not need any modification. The later 

two possibilities are modified to G. . = 1 and G. . = 0 

ij 

respectively. C <’ 0 implies that there are very remote chances 
of tool-failure and C 1 ixaplies the certainty of tool failure 

X J 

before the scheduled hour. 

Thus, the number of failures of a tool, IE. . f is deter- 

13 

mined by the multiplication of G. . and ETS? the number of scheduled 

13 

intervals, SFS and ITE. . can be expressed mathematically as j 

J 

NS = (4.10) 

“ij - r''s-=ij 3 (4.11) 

Thus, the cost expressions are developed on the basis of the above 
assumption. 

Machine and Operating Cost 

This cost component is v/ritten in the sane way as in 
Case 1# Therefore, 

C^ = TG . EJ . (C^ + C^) ( 4 . 12 ) 

Tool Cost 

In case of scheduled tool-replacement cost is incurred 
due to two types of tool-change activity s 

a) tool changing cost becmse of tool-failure before 
the scheduled interval, CGH^, and 

b) the dost becaise of scheduled group replacement, 

CGHg. 
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Tlie for CCH 


^ and OCE^ are ’written as follows ; 


K I. 

c®:^=S S"" 

i=1 j=1 


\ ( Hr. . + "I i 


H n. r- 




i HR. 


ji. 


i=1 0=1 ■*■ ■' 


i Hi 


-ii 


i s 


^ n.. + 1-^ 'i/* "^13 • 


10 


H 


1} 

-J 


EC + .> . ' j IR. . I. TS. . 

10 -0 


. cs 


+ TS . G, 


H ■ H 




IS 


(4.13) 


I 


n. 

X 


i= '' l'=1 L ^13 1 


+ 1 / , EC, , + / 


... \ ' I i isl - 

10 \ L -i 

i^1 0=1 '■ 


IS 


m. 


I n. 

-■. —.1 r- 

; '3 ' r KTsl. 


+ 1 + 1 M . TG. . . CG + / > 1 IS I. TS. . . CS 

13 ^ 1 / 10 4-, , L 13 


7 


i=1 0=1 


(4.14) 


The tool cos.t, OCH, is sum of GCH^ end CGH^. 

Ion-Productive Cost 

Elis cost component occurs for both (a) end (b) type of 
tool-change activify involved and the expression is written keeping 
in mind the same assumptions as in Case- 1. Therefore, 


N n. ^ _ 

’‘nT*""' j' 

CIP = ..7 7. I -=^ 1 . CP . IR +1 ^ . IS . CP 

3^1 3=1 L, 1^ J ' 


i3 i TG J 


■ I 

T 


n. 


. )> TS. . . IR. .■ . C^ 

1 10 10 1 


+ GET . IS . C. 


(4.15) 
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viiere 

GET - time taken, for replacing all tools together, min. 

Total cost can be expressed as the summation of equations 
(4.12), (4,13), (4.14) md (4.15). 

4.3 Constraints 

The total cost of machining, thus deteimined, should be 
optimized with respect to machining .parameters i.e., speed and feed 
rate. However, the choice of these parameters is restricted by the 
number of constraints which result from consideration of the phy- 
sical process of machining and limitation of the machine tool 
under consideration. For exaaple, the spindle speed and feed must 
lie within the range of maximum and minimum speeds and feeds 
available on the machine tool 5 the cutting force daould not ex- 
ceed the maximum force that can be resisted by the work and tool 
holding devices; the cutting po?rer required should be less 'than the 
net power available from the driving motor of the machine tool. 

The machining paraaeters detenaine the tool-life, the 
cutting -forces and the povrer ccnsumed. In most of the literature 
on metal— cutting, the coefficients and the exponents involved in 
tool-life equation aid force equation are ccnsidered to be deter- 
ministic, Iwata (25) has advocated that the deterministic values 
for the constants end Ihe exponents are insufficient for the pur- 
pose of estimating tool-life aad cutting forces experienced by the 
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tool. In the present woik tlie stochastic ciiaractei'istic of the 
exponents and the constants aie taken into account with the help 
of chance constrained piDgramming concept as suggested by Bracken 
and McCormick (6), The concept unriddles the problem of convert- 
ing probabilistic constraints into deterministic ones. In brief, 
the ccncept is as follows ; 

"in optimization problem can be stated as : 

Determine (j = 1, 2, ... , n) to minimize the 

objective function 


P = F (x^, Xg, ... , x^) 

subject to the constraints 
n 

a.. X. b 5 i = 1, 2, ... , m 

3=1 ^ ^ ^ 


vhere 

a. .'s are the technological coefficients and, 

-*-0 

b^'s are the constraint coefficients. 

In this formulation a. .'s and b.'s are random variables mth 

normal distribution, T^hile F is deterministic function of x.'s. 

J 

Such a problem can be solved by the chance constrained progra- 
mming concept. 


Under the chance constraints 


n 


Prob ( T; a. . X. - b ) :^ p , 
j=1 J 1 ^ 


i — 1 5 2 5 




m 


(4.16) 
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find X . ' s to minimize the objective function 

J 


F 


F X 2 , ... ,3!^) 


T/tiere p^'s in equation ( 4 . 16) are predetermined probability- 
levels at vdiich the corresponding constraints must be 
satisfied. 

The chance constraint in equation (4.16) is equivalent to 
the following deterministic constraint as derived in appendix A t 


n 

r: 


n 


m.. X. + Kp .) ( y 

3=1 • 35^=1 


p ^~a. . ^a., X. 

* a. . a. 


- 2 


n 

3,1?:=1 


ij ik 

2 1/2 


13 xfe "3 ^ 


iO 


a. . b. 

13 1 


^a. .'^b. 
13 3 


2 :. + b / ^ b. , i = 1, 2,,. 


( 4 . 17 ) 


■where 


2 


/ 00 o 

J e " dt s 


a. . and b. are the mean values of a. . and b. respectively. 

ly X 13 x- . 

“^a. . and ^b. are the standard deviation of a. . and b. 

13 a- 13 1 

respectively, and 

P and P -u are correlation coefficients 

' a. . a., ■ a. ^ b 

13 xk; X3 X 

between a., and a., and bet-ween a., and b respectively, 

13 XK Xj X 

The chance constrained problems is, thus, reduced to iiie 
equivalent non-linear programming problem to minimize the objective 
function f under the deterministic non-linear constraints in 
equation (4,17).” 
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The constraining factors considered are summarized as 

follows : 

(1) Constraints determined by machine-tool specifications; 

(a) the maximum cutting speed constraiit, 

(b) the minimum cutting speed constraint 5 

(c) the maximum feed constraint, aud 
^d) the minimum feed constraint . 

(2) Constraints determined by machine tool dynamics; 

(e) the allowable maximum cutting-force constraint; and ' 

(f) the allowable maxiinum power consumption. 

(3) Constraints determined by tool— life consideration; 

(g) tool life should be greater than or equal to time required 
to produce a number of jobs. The number is also to be 
determined. 

(4) Constraints determined by the cycle time; 

Ch) the machining time of a cutting-tool shoud be less than 
the cycle time. 

Cutting-Tbrce and Power 

The allowable cutting-force may concern the tool-work 
deflection, ten5)erature generated at tool-work interface and so on. 
However, in the present work, the values of maximum allowable force 
are chosen arbitrarily. The cutting-force and power consumed are 
expressed in terms of cutting parameters as follows : 
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FG 



=<- -^2 ^ 
S f d 


( 4 . 18 ) 


FG « S 
33000 


viiere 


( 4 . 19 ) 


S - cutting speed, fpm, 
f - feed rate, feed/revolution, 
d - depth of cut, 

FG - cutting-force, lbs, 

FG - power, HP, 

Q<o» "^3 noimally distributed paraneters with 

j ^ , 'Xg 5 end as their me sn value and ^ » 

and x ss their standard deviabions respectively. 

Hence, the chance constraints for cutting -force and power become 


prob (fG ^ ? 

prob (PG X PC ) X Po respectively. 

^ ci 


vdnere 

FG - Maximum force allowable, 

max 

PC - Maximum power available, and 

max ^ ’ 

p , p are predeteimined confidence level®, , 

[ c. ■ ■ . 

Pool Life Constraints 

Phe empirical relation for tool - life is taken of the 


fom 





( 4 . 21 ) 
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viiere 

- tool life, min. 

’ '^g.jsnd are consideied as normally distri- 
buted paraneters with, •'■>4^ , , ‘=>*.g jand ^ as their 

meal values and CT” » ? 'T' 9 ad-h rr-, as their 

=^4 ''■'==^5 6 7 

standard deviation respectively. 

The tool-life constraint is imposed becarse it is hi^ly 
desirable that a freshly reground tool last as long as possible. 
However, since tool failure phenomena is inherently probable, the 
constraint can never be satisfied with absolute certainty. Hence, 
the chance constraint can be eirpressed as 

prob ( EP , machining time) p^ 

vihere 

p^ - predetermined confidence level, and 

HP - number of components to be machined before the 
tool failure. 

Cycle Time Constraints 

The constraints regarding cycle time are presented in the 
form as suggested by IPz (22). The constraints take sh^e as 
follows ! 

t.. - TC + 6 ^0 , i =1j2, ... E, 

1-1 

j = 1 , 2 , ... n^^ , (4.22) 

^ is a small positive number and it is essential to 
introduce it to avoid any constraint becoming equal to zero. 
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If ^ is omitted, then at least one constraint will be active end 
the suggested solution algorithm mil stop working. 


The constraints (a, b, c, djandh) are deterministic yiille 
the constraints (e, f 5 and g) are chance constrained and they can be 
converted into deterministic with the help of equation ( 4 . 17 ). Ibr 
simplicity the correlation coefficients in equation ( 4 . 1 ?) are assumed 
to be zero. The original constraints after logarithmic transformation 
of variables are linear, hovrever, the converted constraints are non- 
linear because of the nature of coiversion applied. The resulting 
constraints are expressed as follows : 


Maximum speed constraint - 


g, = ln(Y )“X. . ^0 
^1i ' max'^ 1x V 


Minimum speed constraint 


(4.23) 




2i 


(4.24) 


Maximum feed constraint - 


c = In (f ) - X- . . ^ 0 
’3ii ^ max^"^ 2ii 


Minimum feed constraiit - 


. = -In (f - ) + x„. . ^0 

= 4 ij 213 / 


Maximum force constraint - 


n.. = In (pc ) (lii 10.) 2 : • - 

^3 


(4.25) 


(4.26) 


o< x„.. In (d..)+ ■) isl ■ ■ ) ^ 0 

2 id 2 i 3 313 13 ^"^113 5 x 3 

( 4 . 27 ) 
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® 51 d -of.. 

'^213 


v>^. • . 

, 013 


'lij 


+ 


-< 3j_^ (in d^ ) 


( 4 , 28 ) 


Tool-life Constraint - 


^6i3 




5i3^ ^1i 


^ "^6i3^ ^213 + "^-713 <iij) - IJi (m) 


+ ^ (p ) (s ^ 1/2 


(4.29) 


^ 4 _.. + 0 ^ 


413 

u-j (ind. .)‘ 

"^7i3 


^ 5 i 3 ''^^■"'^"<<613 ""213 


(4,30) 


Cycle Time Constraint — ■ 


®7id “ - ^ (^iP + ^11 + ^213 + Ui (TO) + e ^0 (4.31) 

Humber of Components Constraint - 


gg = “ In (m) + 3n (IJ) X 0 


(4.32) 


viiere 


x^. s= In (y. ) , X = In Cf 
■ji Vi^ » ■2i3 


The above constraints are for all i = 1 2 


H and 


3 3 y 2 y • , , y n^ • 


In the above constraint set the index i refers to spindle 


and the index 3 is for an operation on i"*^ spindle. It is to 


be 
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noticed that the above constraint set is written in a general foim 
applicable to both type of automatics under consideration. 

The constraint set presented above does not include the 
power constraint. It is becaj.se the stiucture of this constraint 
will be different for two types of automatics under consideration. 
The horse povrer constraint for the automatic type - A is as follows; 


Sgi 


^i 

In (PC )- 'V'lnPC + V-^Cp^.) (g’ > 0 


n. 


g'q, = y (PG^ i In (10.)^ , . (T~ 


i=i 


13 


c< 


1i3 


(4.32) 
2 


1i 


2i3 


+ X, 


■2ij 


cn^ + (in d- ^ 

3ij 4.ij 


(4.33) 


for all i = 1, 2, , II. 

Hence, there will be as many constraints of type (4.32) 
number of spindles. 


The horse power constraint for the automatic type - B s 

In this case there will be only one constraint. 

1 n. 

(4.34) 


H n. ■ 

n •i -1— '4'' ■ " V ' ’ 


g 


i=1 3=1 


?r 

013 


+ ^ . .til + (itt d. .) 

11 213 


(4.35) 



46 


The total number of constraints, for example, for a 
case of five spindles automatic type - A and n^ vector (2, 1, 1, •Jj)j 

vail be 44, vhile, for type - B the number is 30. The variables 
correspond to cutting parameters. This number will be 11 and 4 in 
tvro types of automatics respectively for the above example. 

Thus, one finds that the problem reduces to an optimiza- 
tion problem with non - linear objective function constrained by 
non-linear constraints. The available mathematical programming 
teciiniques ensure of a global minima only for the conVex functions. 

But both the objective flinctions under considerations are not 
convex md none of them is having a expression for partial deriva- 
tives, therefore, any optimum solution obtained can not be claimed 
as global minima. 

Chapter V discusses the solution algorithm of penalty 


function ^pioach. 



CHAPTEE Y 


OETBIIZAEIOI PROCEDUEE 

5.1 Int ID duct ion 

The optimization piobleia formulated in Chapter lY in its 
generalized form can be stated as ; 

Minimize f (x) 

s»t, (x) ^ 0 j i = 1j 2} 

■where x is the vector of indep ardent variables (x^, x_, x ) 

corresponding to cutting parameters. The constraints form a bounded 
region in n - dimmasions, 

A number of non-linear prograEaning methods such as the gra- 
dient projection method (40), Zoutendizk’ s feasible direction method 
(5l)j penalty function method (is), etc. have been discussed in the 
literature. Most of the suggested methods axe of fairly recent origin 
(l958 and onwards) and a fool proof answer to the problem of selection 
of a method is not yet establi^ed. Ho'wever, tbe experiences of 
different, research workers is always of help. Pox (22) reports 
that the gradient projection method has been described for the gene- 
ral non-linear programming problem, but its effectiveness is mainly 
limited to problems in which tie constraints g^ (x) are linea 3 >- 
functions of x . lasdon (3l) opines that interior penalty function 
approach is particularly attractive in dealing with problaas that 
have markedly non-linear constraints and the objective function 
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since, it approaches the solution value from inside the constraint 
seo, and thus, avoids the difficulty of movement along non-linear 
bomdaries of the set. -vihich is required in the other methods* 

Hence, the interior penalty function method is adopted in this study. 

5.2 Optimization Methodology 

The optimization procedure consists of converting the 
general constrained problem to a sequence of unconstrained problems. 
The sequence of such problems are solved by Davidon - Pletcher - 
Povrell's Variable Metric Method (2l). Golden Search is employed to 
carry out cne variable minimization. 

5.2.1 Conversion to Unconstrained Minimization 

Piacco and McCormick (l8) have developed an algorithm for 
transforming the general ccnstrained problem involving inequality 
constraints to a sequence of imconst rained minimization problems. 

The transformation defines anew function ( x , r ), 

m 

0 ( z , r ) = f ( X ) + r . ^ ^ ( 5 . 1 ) 

and the minimization is carried over a strictly moriotonic decrea- 
sing sequence of r - values. 

Piacco and McCormic have shown that if 

(i) the interior of the constraint set is non empty, 

(ii) the function f ( x ) and g^ ( z ), (i = 1, 2, ..., m) should 
have second derivatives continuous, 
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(iii) the set of polats in the constraint set for which f (x) < Y 

^ 0 

is bounded for every finite T , aid 

o 

(iv) the function f ( x ) is bounded below for x in the -cons- 
traint set, 

thenj; the optimal solution to the unconstrained problem 
app3X(aches a local minima of the constrained problem as the 
value of r approaches zero . If ,, in addition 

(v) f ( X ) and ( X ), (i = 1, 2, m) are convex functions, 

and 

(■vi) ( X, r ) is strictly convex in the interior of the cons- 

traint set for every r > 0 , 

■thoi, the optimal solution to the unconstrained problem 
^preaches the absolute minimum of the constrained problem as 
r approaches zero. 

The steps comprize of ; , 

1. Start with x , an initial solution vector, which 

0 

must be strictly inside the constraint set, and 
r^ 0* let i~'l,2,**«,n* 

2. Minimize 9" ( > ^i ^ starting from x^ ^ , aid 

subject to no constraints^ 

5. If the convergence criteria is satisfied, go to 

Step 4, otherwise ■reduce r by choosing r^^^ < r^ 

aid return to Step 2 with i replaced by i -t- 1 . 
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4. T^e r = 0 and minimize the objective function and the 
iteration is teuninated when the convergence criteria 
are satisfied in -the unconstrained minimization too tine. 

There are number of points to be considered in applying 
the method : 

a) The starting solution required by Step 1 has to be 
Carefully chosen, i,e., not very far from optima to 
avoid excessive computational efforts, 

b) Care has to be exercised in the initial value for r, 

c) The convergence criteria in Step 3 has to be carefully 
selected. 

In tie present piobleni, x^ has been chosen by a carefUl 
judgement. The convergence criteria chosen is simple. The Step 3 
is terminated vdien ihe difference in tvro successive finctional values 
is relatively small* This is a weak criterion - a better one being 
a check on the difference between the values of x^_^^ and x^* However, 
in the present usork the first criteria is considered. 

Ibr ihe selection of initial value of r, a method sugges- 
ted ty piacoo and McCormick is explained in Appendix B, 

5*2,2 Procedure for Unconstrained Minimization 

An unconstrained minimization of the function <5^(x, r) 
must be carried out over Ihe design vector x Ibr each value of r 
used in Hacoo - McCormic algorithm. The melhod makes use of Ihe 



quadratically convergent iterative descent: method stiuctured by 
Davxdon aid subsequently modified by Pletdier and Powell (2i). 
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The method uses the concept of local Hessian j"'' , but 
approximates it by a metric The procedure for computing the 
metric completely eliminates the need for evaluating second deri- 
vatives and performing matrix inversions, and yet the metric vdiich 

is improved at each iteration converges to j""* . 

m 

The steps in the algorithm are : 


1. Start vdth a positive definite matrix (usually ciiosen 

as the identity matrix) and an initial point x . . The 

.th , . ° 

1 step, 1 = 0,, 1, 2 ... proceeds as follows : 

2. Compute the gradient vector '\7'r’(x,r) = G 

i i 

3. Compute the direction "s^ = ~ x. , r ) 

. _ * 

4* Ciioose a step length ‘X ^ to minimize 



P ( ) = 95(x + o<'. S., r ) 

/ ' 1 'll’ ’ 


Compute ^ i 


Compute anew value x. from the relationship 


i+l 

x.^. = X. + C< . ’s. 

1+1 1 1 1 


Compute ( X. 


i+1 ’ 


) = G. 


i+ 1 


Compute y^ = - 


Compute the matrix 


C£NTfrA). -f/ '[ 


H 


Sm, J 


26289 
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10* Compute the matrix B. 


y'i H. y. 


11# Compute the successor 



H, + A. + B. 
1 X 1 


12, Check for the convergence criteria, if it is satisfied, 
return to constrained minimization routine. Otherwisie, 
return to Step 2 using the successor matrix as the new 
and replace i by i + 1 . 


In the above algorithm the gradient vector (calculated in 
Steps 2 and 7) is found out. The method of finite differences is 
employed here and is discussed in the Appendix C. The Step 4 is ex- 
plained in section 4.2.3 of this chapter. 


In applying the method care has been exercised to ensure that 
the H - matrix ^ould always be positive definite gad it is not 
updated with data arising from poor approximation to Ibx has 

mentioned many ^pioaches to this problem. Jirst, the algorithm 

* T 

used for computing . msy be reapplied until S: G. is suffi- 
ciently small, another alternative is to simply skip the updating 

m 

Steps 9, 10 and 11 of the algorithm vfaen S G. . is too large. In 

1+1 

. .. ■ ■ ■ ■ . ■ ■ , , . ■ ■ ■ 

other words if o< . is not close enough to minimum along S. , H = H 
1 1 ’ i+ 1 i 

and = ^i +1 *^i+-i process is continued. 

It is difficult to choose between these approaches. The 
first metiod may require excessive computations to refinecs< 
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at points far from but the second approach may skip valuable 

opportunities to improve the H matrix* A compromise is struck at 

this state; 'X is calculated twice and if value of G is 

^ 1 i+1 

still la^rge^ the updating steps are skipped* In the present work 
another method suggested by Pox is employed. A number C = G / 

m 1 x+1 ^ 

,T 

calculated. The contention is'liiat C = 0 at points of 
mininMn. At all points along in the neighbourhood ofcK*^ , G«l1, 

The value of 0 chosen is .001. 'So if the calculated C ;?-,001, tiie upda- 
ting steps are skipped. 

Another area of numerical difficulty is the classical round 

off problem. If the scsling is bad| depending on the scale factor, 

or = A^ md the latter form is singular. The 

remedies are : increase precision of arithmatic ora proper scaling 

of the variables, or continuous check on H matrix. The last method 

is put to use in the progranme. The remedy consists of setting H 

matrix to or some other predetermined positive definite matrix 

and the process is continued. Thus a check is continuously kept on 
T 

the factor G^. iny time it is positive, the matrix H is 

T 

reassigned the value of identity matrix. If G^ G^ is at positive 
level, it implies that the direction - G^ is no longer a direc- 
tion of descent and H matrix has to be modified. 

In the algorithm are normalised to avoid excessive 
increment in design variables at any iteration. It is possible 
becaise in the present problm,1he variable are transformed by 
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taking logarithm viiile is at nomal sc^e. 

consists of = s. / (s.s)'*'^^. 


The normalisation 


The termination criterias consists of check on step size 
end 3 ^. If g^y ^f ^ msgnitude less than a pre- 

a^signed Cr , ( ^ > o, e«i ), the process is terminated. 


5.2,3 Golden Search iOgorithm (18) 


The numerical difficulty ejcperienced with inaccurate o(^ 

i 

hao alreadj’’ been discussed. The field of uni-directional minimiza- 
tion with all the sophistication attained remains a comparatiTrely 
crude area. Golden Search method lAhidi is used in the present work 
is one of the best methods available but still relies on arbitrary 
but so called Pibonacci number (l8). 

The steps are as follows : 

1. Pirst an upper bound is obtained. The first lower 

bound = 0. ©^ is obtained bj?' evaluating the 

ftinction (x, r) at successive points whose © values 

are in the Pibonac_ci ratio , 1.618 = (l + /5 ) / 2| 

J 

that is © is ( 1.6 16)'’ 

3=0 

where J is the smallest non negative integer 3 such 

that + J' (l.6l8)^ S. I > ^(x.). 

k=0 "--j ^ 
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2. The interval boundary e is reduced by computing two 

values of © in ihe interval ( © - © ) , 

u L 


© 

a 


h * ( ®u ■ h* 

»L * •«« ( ®u - ®i) 


3* The Ci values of the (interior) correspond to © abd 9 

a t 

are compared* 

4 . If 0 ( + ©^ 5 r) <. + ©^ , r) 5 then 

® ^ ®b ’ of the pioperty that .382/ 

0.618 = 0.618, by letting © = ©,,'©' = © and 

u D ^ b a 

©1 = © and recomputing ©' = © ' + .382 ( ©' - ©’ 

b a I 'u L' 

Step 3 can now be repeated. 


5. If (x^ + ©^ S^, r) > <56 (x^ + ©^ S^, r), then by 

letting ©I = © , ©' = © , end 9’ = ©. and computing 

•U a a c u u 

^b ~ ^ Step 3 cm be repeated. 

6. If (x^ + ©^ S^, r) = + ©^ S^, r), let 

©£ = 0^5 ©^ = ©^ and repeat Step 2. 

7. 75iea ©^ - ©^^ is acceptably small, © is approximated 
by 9 = © + © / 2 and X. . = x. + © S. . 

Is U 1-fl 1 X 

The procedure fbr obtaining -^e optimum -value as outlined 
above ensures a global optimum cn^y for convex functions, Ibr . 
functions not satisfying the critefion the minimum obtained may be 
but a local minimum only (31 ). 



CEAPTER VI 


COlffUTEEIZjfflIOF Mil) TESTffiG OE THE IIODEL 
6.1 Co lapu ter Package 

The optimization procedure explsined in Chapter V has 
been computerized so that it may be applied to the practical 
situation of machining of parts on multi - spindle automatics. 

The program also generates the constraint set (equations 4,23 to 
4.35). The generation of constraint subroutine is quite general. 
The objective functionsj 1) minimization of cycle time and 
2; minimization of total cost of pxoductionj have been progranmed 
The computer listing is put in appendix E. The comment cards are 
pub liberally to explain the ■'.".orking of the program and the mean- 
ing of Various important variables used in the program. 

In tlie case of multi-spindle automatics the amount of 
information to be fed to the Computer is enormous compared to the 
single-tool case. To a great extent this difficulty is overcome 
by representing the information in a matrix form. This method 
is advantageous also in many waj''s ; 1) r’evelopment of a genera- 
lized computer package is possible, 2) use of the package becomes 
easy, 3) the development of logic is simplified, etc.. 

In this work, an. attempt has been made to develop a 
general purpose package. The optimum machining parameters can 
be determined for any operation in Tshich metal is removed in 
form of chips. Eor example, single tool operation on lathe, 
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turret lathe operation, mlling operation, etc. cai be handled, 
however, the package has been tested only for the multi-spindle 
aitomatic operation. The package can haidle m automatic • ■ 
liavii..t3 upto 8 spindles. A aaKimum of 6 operations can be perfor- 
med on each spindle. The machine tool specification determines ijhe 
size of the 30b vtiich the machine cm handle. 

ihe package has been tested for a saaple 30b shown in 
4. The input infonnation needed for the determination of 
optimurii cutting parameters is given belov/. 

6*2 Operation Plan of the Job siiovsn in Pig. 4. 

Spindle 1 - Stepped tuming to produce surface 

( 2 ) and ( 3 ). 

Spindle 2 - Pozm turning to produce ( 4 ) and (s) surfaces. 

Spindle 3 - Porm turning to produce surface (6) and 
spot drilling. 

Spindle 4 - Prilling of hole 1/4". 

Spindle 5 ~ Parting off md thread cutting. 

6,3 The Input 

The infoimation regarding, a metal removing activity 
can be grouped under the following headings ; 

1. Eaw mateiial characteristics. 


2. Component requirements. 
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3. Data about ttie component to be madilneds 

4. Machine characteristics, 

5. Tool characteristics, 

6. Meai and standard deviation of tocJlr-life and 
cutting-force equation exponents, and 

7. Cost and time parsmeters. 

In the follOTong sections the various elanents pertaining 
to a data group along vdth their representative values taken in 
the present wo-ii are described. 

6.5'»1 Raw Material Characteristics 

This group of information is named as Raw Material 
Characteristics Matrix derioted by BlifflR. lable 1 shows the 
RORTRiR name and the element along wilii value . 


TABIE 1 

Raw Material Characteristic Matrix 


RORTRAR 

RJilE 

Description 

Yalue 

MMB. (l) 

Type of material used (Code Ifo.) 

1060 

iMATR (2) 

Brinell hardness of tiie raw stock 

250 

ilATR (3) 

Yield strength of the raw material, psi 

127500 

MATH (4) 

Any important dimension (length, dia- 
meter, etc.) of the material, inches. 

1.2 

AMATR (5) 

Raw material cost rate, Rs. per/unit 

1.00 



6.3.2 Component Requirements 


Tile component requiiements are stored in Component 


Requirement Matrxx denoted as COIffi. The various elements are 
put in Table 2 along with their assumed values. 


TIELS 2 

Component Requirement Matrix 


DORTRj4I7 

imsE 

Description 

7alue ■ 

COMR 

(1) 

Humber of Components to be manu- 
factured. 

10000 

COMR 

(2) 

Material grade specification (Code No.) 

1060 

COMR 

(3) 

Raw Material length required for 
component, inch. 

2.6 

COMR 

(4) 

Attrition percentage. 

5.0 

COMR 

( 5 ) 

Date of Receiving the order. 


COMR 

(6) 

Due date 



6.3.3 Data about the Component 

This information is put under the matrix. Component 
Data Matrix, designated by C0I£D. This is a three dimensional 
matrix. The first dimension is allocated to the spindle, the 
second is for operations, and the third is for various elements. 

The specific values have been generated from the job 
(Dig, 4 ) under consideration. The description of the input 
elements is as follows s 
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The subscript I refers to spindle and J to an operation 

t;ii 

on I spindle. 

1. GOI-ID (ij J, l) - Tool travel length, inches. 

2. GOIID (l, j, 2) - Depth of cut, inches. This is specified only in 

tile case of multi - spindle operation ^^here in 
the job is finidaed in one cycle of the madiine, 
O'tlierv/ise, it is a machining parameter to be 
determined, inches. 

3. GOMD (l, J, 3 ) - Feed indez. It is a logical variable -viiich 

governs the construction of feed constraint. 

For example, if the two tools are employed cn 
the same slide, the feed constraints are to be 
formed only for any one tool. Its value is 1 
if new feed constraints are to be formed* Other- 
wise it is 2. 

4. GOIvID (l, J, 4 ) - Speed index. It is also a logical variable 

governing the speed constraints. For example, 
if the tfjo tools are put on the same spindle, 
the speed in rpm for both the operations has to 
be the same. The constraint is generated only 
for one tool. It is 1 if new constraint set 
is to be generated. Otherwise it is 2. 

5. GOtm (l, J, 5 ) - Machining length, inches. This value indicates 

the actual machined length by a tool. In some 
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cases it may be different than tool travel 
length (COMD (l, J, i)), inches. 

S, GOilD (l, J5 6) - Index for secondary operation or primary opera- 
tions. This is a binary variable. It is 0 for 
a secondars' operation and 1 for a primary 
operation. 

7. COLiD (I5 Jj 7) - Depth of cut index. This index helps in deci- 

ding about TAiaether the depth of cut is a 
parameter to be determined or it is a given 
parameter. It’s value is 0 and 1 for the two 
cases respectively. 

8. ooim (l, J, 8) - It is an index to indicate whether the operation 

is to be done on the axial slide or cross slide 
of the automatic. This binary number also helps 
in constructing the constraint set. Ibr example, 
if on a particular station the two operations are 
being done, one operation on axial slide and the 
other on cross slide, "tiie feed rate can be diffe 
riht for axial slide than for the cross slide. 

It is 0 if the operation is done on axial slide, 
otherwjise it is 1, 



table 3 

COIvIPOlTMT data MAIEIZ 
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(COMD (l. 

J, 

IC), K 

= 1 

5 8) 



Spindle j[Opn. Ho 

4 1 i 

2 

1 3 

I 4 ! 

5 1 

6 

1 7 

1 

8 

Ho, 



\ inch 1 

Jnch 


1_ 

1 inch . 1 


1 

JL 



J = 1 

1,250 

0.125 

1 

1 


1.250 

1 

1 


1 

I = 1 













J = 2 

1.250 

0.125 

2 

2 


.570 

1 

1 


1 

1 = 2 

J = 1 

0.500 

0.870 

1 

1 


.500 

1 

1 


1 


J = 1 







0 




1 = 3 













J = 2 

0,125 

0.250 

1 

1 


.125 

1 

1 


1 

M 

11 

4^ 

J = 1 

2.500 

.125 

1 

1 


2.500 

1 

1 


0 


J = 1 

0.125 

.440 

1 

1 


.125 

1 

1 


1 

1=5 













J = 2 

— 

- 

- 

- 


- 

0 

- 


- 


6,3»4 Machine Characteristics 

This group of information is fed to the computer under the 
matrix name Machine Characteristic Matrix (mCSI) . It is a two dimen- 
sional array. The first subscript relates to the number of madiine 
tinder consideration and the second subscript is for various elements. 
The table 4 shows the description of various elements. Here, the 
characteristics of one machine are ^own, because only one machine 
is taken and the sensitivity analysis is carried out for its horse 
power caqpacity only. 
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TiffilE 4 

MACHINE CHARAOTEHISTIC MAI3IZ 


EOEETiAI 
1=1 

SME 

i Description 

Value 

aiG (i, 

l) 

lumber of stations on the machine 

5 

CMC (l, 

2) 

lumber of axial slides 

5 

CMC (l, 

3) 

Type of madiine (Tjrpe - A or Type - B) 

(1,0) 

CMC (l. 

4) 

Maxiaum speed available, rpm 

2400 

CMC (l, 

5) 

Minimum speed available, rpm 

100 

CMC (l, 

6) 

Minimum feed rate, ipr. 

.0005 

CMC (l, 

7) 

Maximum feed rate, ipr. 

.2000 

CMC (l, 

s) 

Minimum depth of cut, inch 

.0100 

CMC (l, 

9) 

Maximum depth of cut, inch 

.5000 

CMC (l. 

10) 

Horse power 

20.0000 

CMC (l, 

1l) 

Power Transmission efficiency 

.8000 

CMC (l, 

12) 

Axial minimum feed rate, inch/rev. 

.0005 

CMC (l, 

13) 

A>:ial maximum feed rate, inch/rev. 

,2000 

6.3.5 Tool 

Characteristics 



A cutting tool, generally, associates itself with 
a lazge number of specifications. In the present work only the rele- 
vant characteristics are specified. Ihis group of information is 
put under Tool Characteristic Matrix denoted as TCSI, The descrip- 
tion of elements is given as follows. 
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1 . 

2 . 

3. 

4. 
5 » 
6 . 


TCIvI (l, J, l) - 
TCM (l, J, 2) - 
TCM (l, J, 3) - 

TGM (l, J, 4) - 
TCS'I (l, J, 5) - 
TCII (is J, 6) - 


lumber of teeth on the tool. 

Diameter of the tool (if releTan.t)s inc±i. 

Maximum force or torque allov/able lbs. or lbs. in 
respectively. 

lumber of regrinds permissible 
Tool edge cost, Rs, , and 
Horse power constant. 


The values corresponding to the elements of matrix 
TCM is shovai in Table 5. 


TABLE 5 


Tool Characteristics Matrix 



“T 

_ i 



Tdf (l, J, 

K)s K = 

Is 6 


Spindle 

I f 

V Opn. lo.v 

1 

1 2 


4 I 

5' 1 

6 e 

Ho. 

_L_^_ A 


f in. 

I lbs. 1 


Rs. i 

X 10^ 


J = 1 

1 


200 

3 

2 

.795 

I = 1 

eg 

II 

1 

- 

200 

7 

J 

2 

.795 

1=2 

J = 1 

1 

- 

200 

3 

10 

.795 

11 

. 

J = 1 

- 

- 

- 

- 

- 

- 


J = 2 

- 


100 

3 

5 

.795 

1 = 4 

J = 1 

2 

• 25 

50 

3 

3 

7.950 


J = 1 

1 


50 

3 

5 

.795 


1 = 5 
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6.^.6 Exponents of Tool-life and Cutting-Eorce Equations 

Tile mean values of tlie exponents of tool-life and cutting- 
force equation for tlie cutting tools involved in the example are 

taken from Kroiienburg- (30). The standard deviations of the exponaits 
are chosen arbitrarily. 

The meal values and standard deviation of exponents for tool- 
life and cu btin^ force equation are tiut in Table ^ and 6 respectively. 


TABLE 6 

Cutting-Eorce Equation Exponents 


Spindle fOpsratior 
Ho. 1 Ho. 

1 _ 


MEM fALUES 


'1 

stJdaed ieviatioi 

TMateriai 

i-loj 


oT ' 
1 -=<2 1 

'^3 1 

"k 0 

Tr~r 

'"“■2 

T 

1 


1 

!S 


j = 1 

3.9 

0.0 

0.8 

1.2 


0,4 

0.0 

0.2 


0.2 

JL 

HSS 

I = 1 














j = 2 

3.9 

0.0 

0.8 

1.2 


0.4 

0.0 

0.2 


0.2 

HSS 

T ~ 0 

-»■ > 


U.O 

0.8 

0.0 


0.4 

0.0 

0.2 


0.2 

HSS 


J = 1 


mm 










1 = 3 














J = 2 

3.9 

0.0 

0.8 

0,0 


0.4 

• 0.0 

0.2 


0.2 

HSS 

1 = 4 

J = 1 

3.9 

0.0 

0.9 

0.0 


0.8 

0.0 

0.2 


0.2 

HSS 

I = 5 

d = 1 

3.9 

0.0 

0.8 

0.0 


0.4 

0.0 

0.2 


0.2 

HSS 


CM 

il 

- 

- 

- 

— ^ 


— , 


' M ' 





^ •3ml Time and Cost Parameters 

The representative values of cost parameters are as follows : 
C_| - Direct labour cost rate = .02 Rs./min. 

Cg - Oveihead cost rate . = ,01 R3./min. 
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=3 

- Machine setting cost rate = 

.05 Rs./min 

GG 

- Grinding Id^our cost rate = 

.03 Rs./min 

C£ 

- Tool-change cost rate = 

.05 Rs./min 


Time parameters are different for different tools. The ele- 
ments are grouped under matrix TEIEG, The description of the elements: 
is as follows : 

TBIEC ^(ij J} l) - Tool - Change time, min., and 
TBIEG (l, J, 2 ) - Tool - grinding time, min. 

TABLE 7 

Tool-Life Equation Exponents 


SpindT^6p¥ration"jr _ mm YALUES 


STMLAED EEYIATIOIT 


Ho 

1 » 

JL 

Ho . 

5 41 

^1 

5i 

6 1 

X 

7 ! 

‘^4 5 


a~:“r 
' 6 J 





J = 1 

8.9 

CM 

'=3- 

1 

-1.8 

- 0.5 

0.25 

0.20 

0.05 

0.03 

I = 

1 














J = 2 

8.9 

-4.2 

-1.8 

- 0.5 

0.25 

0;20 

0.05 

0.03 

I = 

2 


J = 1 

9.5 

-4.2 

-1.8 

- 0.5 

0.25 

0.20 

0.05 

0.03 




ii 







_ 


I = 

3 














J = 2 

9.9 

-4.2 

-1.8 

- 0.5 

0.25 

0.20 

0.05 

0.03 

I = 

4 


J = 1 

10.3 

-4.2 

-1.8 

- 0.5 

0.25 

0.20 

0.05 

0.03 




J = 1 

9.9 

-4.2 

-1.8 

- 0.5 

0.25 

0.20 

0.05 

0.03 

I = 

5 ' 














J = 2 










The table 8 shows the Talues of -the elements in matrix 


TUfflC 
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TiBEE 8 
Time Parameters 


Spindle Bo. 

1 Operation Bo. 1 TIFJIC (l, J, l) J TDIEC (l, J, 2) 

_i i Min. f Min. 


J = 

1 

20.0 

5.0 

1=1 


J = 

2 

20.0 

5.0 

1 = 2 

J = 

1 

25.0 

10.0 


J = 

1 



1 = 3 


J = 

2 

25.0 

10.0 

. 

1 = 4" 

J = 

1 

15.0 

20.0 


J = 

1 

10.0 

10.0 

1=5 


J = 

2 

- 

- 


The above set of data is fed to the computer as part of 
the main program. The program geners'tes the constraint set and 
evaluates, the objective function. The initialization of all the 
cutting parameters is done in the constraint generation program. Sub^- 
quently the logic enters into the SHE algorithm aid finds out 
the optimum solution. 



CHiPTEE 7II 


EBSUIUS MD OOIGinSIOHS 

7.1 Eie Ee suits and Conclusions 

The results obtained for the example (Pig. 4) des- 
cribed in Chapter VI are discussed here. 

The optimun cutttog conditions and. the corresponding 
values for the objective functions are obtained. The optimiza- 
tion is carried out with respect to the minimization of cycle 
time and minimization of total production cost. The results 
are tabulated in Tables 9 - 12. The tables 9 and 10 corres- 
pond to the minimization of cycle time for a machine - tool ha- 
ving edacity of 20 and 10 horse power respectively. The tables 
11 and 12 correspond to ihe production cost minimization aid 
show the machining paraneters and total cost incurred for un- 
scheduled and scheduled tool-replacement (2-A) procedures res- 
pectively. 

Both type of automatics - the automatic having diffe- 
rent speed aid feed on each spindle (Type-A) and the automatic 
having common speed and feed for each spindle (Tybe-B) -have 
been studied. The results for both of them are tabulated in 
Tables 9-12. 

The nomenclature used for the rows of tables 9-12 
is explained in Table 8a* ^ 

The efforts are directed to study the effect of diffe- 
rent probability levels concenaing the constraints (equations 



TiBLE 8 a 


x11 
x12 
x13 
x14 
x15 
x21 
x22 
x23 
x24 
x25 
. HP 

til 

tl2 

t21 

t3l 

t41 

t51 

Off 

Til 

T12 

T21 


Table Of Nomenclature Eor Tables 9} 10, 11, and 12 

Eeed rate for first spindlej ipr. 

Feed rate for second spindlej ipr. 

Feed rate for third spindle, ipr. 

Feed rate for fourth spindle, ipr. 

Feed rate for fifth spindle, ipr. 

Speed for first spindle, rpm. 

Speed for second spindle, ipm. 

Speed for third spindle, ipm. 

Speed for fourtli spindle, rpm, 

■.S^«ad for fifth spindle, ipm. 

Humber of components for -which -the tool-life constraint is 
to be satisfied. 

Travel-time for tool 1 on spindle 1, minutes. 

Travel time for tool 2 on spindle 1, minutes. 

Travel time for tool 1 on spindle 2, minutes. 

Travel time for tool 1 on spindle 3, minutes. 

Travel time for tool 1 on spindle 4, minutes.. 

Travel time for tool 1 on spindle 5, minutes. 

Cycle time, minutes 
Iiife of tool 1 on spindle 1, minutes, 

life of tool 2 on spindle 1, minutes, 

life of tool 1 on spindle 2, minutes. 



2)31 

T41 

T51 - 

NS 

SL 

PT 

PC 

TOC - 
TCC1 - 

rcc2 - 

ITC - 
ITC1 - 
I!rC2 
TO 


Table 8a (Contiaued) 

Life of tool 1 on spindle 5, minutes. 

life of tool 1 on spindle 4} minutes. 

life of tool 1 on spindle 5j minutes. 

lumber of tool-replacement schedules in scheduled tool- 

replacement procedure. 

life of tile system, unscheduled tool-replacement procedure, 
minutes. 

Total machining time, minutes. 

Cost corresponding to machining time, Rs, 

Tool cost in case of unscheduled tool-replacement procedure, Rs. 
Tool cost in case of scheduled tool-replacaaent for the tools 
changed because of tool failure before the scheduled inteiival* 
Tool change cost due to group replacement in scheduled tool- 
replacement procedure, Rs. 

Idle time cost corresponding to TCG| Rs. 

Idle time cost corresponding to TCG1, Rs. 

Idle time cost corresponding to TGC2, Rs. 

Total cost of production, Rs. 



TJffilE 9 
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Cycle Time Itflinimization - 
Horse Power 20 


. ? . | p£:A H ORSE' PQ¥EB=20 ~ jgJTOMAT IC ‘ TYPE'-B HOHSE "PO’ffilfcrPO . 0 




=.5000 


3x = p 


.8413}=. 


P^1 ^"^2^ 

\ - 4 \^ 


PoilP. 

^ 3 }’ 


!2?1 


= .9987 .50OO1 


^ 3 ! 

=.84l3^y= 


'1 “ ^2 
= .9772^ 


- ^3 

.9987 


1. xlixio^ 

2, x12x10^ 

•z 

7.9 

19.7 

14.1 

19.1 

5.6 

7.3 

11.2 

12.1 

4.3 

3.4 

2.7 

4.3 

5.6 

7.3 

11.2 

12.1 

3. xlSxlO-^ 

Z 

1.3 

1.2 

17.9 

3.4 

5.6 

7.3 

11.2 

12.1 

4. xl4x10'^ 

Z 

7.8 

15.8 

8.3 

20.2 

15.2 

16.2 

23.7 

2J.1 

5. .X15X10-' 

6. x21x10"'' 

7. x22x10”'’ 

8. x23x10”'^ 

9. x24x10 

10. xPSxlo""’ 

11. .iTOxlo""^ 

2.1 

2.1 

1.4 

3.9 

5.6 

7.3 

11.2 

12.1 

78,0 

56.9 

81.7 

145.7 

108.9 

145.3, 

239.9 

239 .. 9 

6 1.8 

85.0 

112,4 

167.7 

108.9 

145.3 

239.9 

239.9 

92.9 

142.9 

105.1 

134.5 

108.9 

145.3 

239.9 

239.9 

126.9 

119.1 

178.8 

239.9 , 

108.9 

145.3 

239.9 

239.9 

79.1 

124.2 

156.0 

16 1.6 

108.9 

145.3 

239.9 

239.9 

20.5 

37.4 

47.0 

25.1 

20.0 

20.0 

20.0 

20.0 

12. tllxlO^ 

2.0 

1.1 

1.1 

0.4 

2.0 

1.2 

0.5 

0.5 

13. t12x10^ 

2.0 

1.1 

1,1 

0.4 

2.0 

1.2 

0.5 

0.5 

14. 112 1x10 

1.3 

1.2 

1.2 

0.5 

0.6 

0.3. 

0.1 

0.1 

15. t31x10’’ 

1.7 

1.1 

0.1 

0.4 

0,3 

0.2 

0.1 

0. 1 

16, Mlxio'^ 

2,5 

1*5 

1.7 

0,5 

1.5 

1.1 

0.5 

0*5 

17. tSlxlo"^ 

1.8 

1.2 

1.4 

0.5 

0.5 

0.3 

0.1 

0.1 

18. CSJxlo'' 

2.5 

1.3 

1.7 

0.5 

2.0 

1.2 

0.5 

O’. 5 

19. Tllxio’’ 

62.8 

84.8 

55.2 

4.6 

22.9 

7.6 

1.2 

1.1 

20, 

21. T2lxic“^ 

c, o n 

S^ A n 

■ZVi 0 

/ c 

^0 Q 

a 

. 0 

-1-4 

93.3 

41,2 

20.3 

4.1 

37.0 

12.4 

1.9 

1.8 

22. T31x10“^ 

143.1 

28.3 

9,6 

15.5 

108.1 

36.1 

5.5 

5.4 

23 . 141x10""^ 

51.2 

36.7 

16.7 

3.4 

20.0 

7.7 

1*3 

1,2 

24. TSlxlo”"* 

187.0 

46.9 

30,8 

12.9 

111.3 

37.1 

5.6 

5.5 
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TJfflIE 10 


Cycle Time ISLnimization - 
Horse Po\’'/er 10 




. TYPE- 

a ': 

HOESE PC 

1 ^ 

= *=3 
=.9772| 

rfEB=10.0TSri0I:;ATIC TTPE-B horse PO\7EE=io.O 

ilP- = P^iP 

I =.5000|= 

1 = P2 

= ^3 
.8413 

¥ 

5 

i. 

Pi = P2 
= P3 

=.9987 

Ip- = P, 

1 “'’3 

1 =.5000 

r 

1 

P- = P 
= P 

=.8413 

1 

’1 = ^ 

= P3I 

=.97721 

P.j = P2 
= P3 

=.9987 

3 

1. xllxlO 

7.9 

19.7 


14.1 

19-1 

5.6 


7.3 


11.2 

12.1 

2. x123c10^ 

4.2 

3.4 


2.7 

4.3 

5,6 


7.3 


11.2 

12.1 

3. xl3x10^ 

*7 

1.3 

1.2 


1.8 

3.4 

5.6 


7.3 


11.2 

12.1 

4. x14x10 

7.8 

15.8 


8.3 

20.2 

15.2 


16.2 


21.7 

22.1 

5. x15x10^ 

2.1 

2.1 


1.4 

3.9 

5.6 


7.3 


11.2 

12.1 

6. x21x10’''^ 

78.4 

56.9 


81.7 

145.7 

108.9 


145.3 


239.9 

239.9 

7. x22x10 ** 

61.0 

85.0 


112.4 

167.7 

108.9 


145.3 


.239.9 

239,9 

8. x23x10"’' 

83.9 

142.9 


105.1 

134.5 

108.9 


145.3 


239.9 

239,9 

9. x24x10"'’ 

128.2 

119.1 


178.8 

239.9 

108.9 


145,3 


239.9 

239.9 

10, X25X10'''' 

80.3 

124.2 


156.1 

161.6 

108.9 


145.3 


239.9 

239.9 

11. HPxlo"”' 

19.8 

37.4 


47.1 

25.1 

20.0 


20.0 


20.0 

20,0 

12. tllxio’' 

2.0 

1.1 


1.1 

0.4 

2.0 


1.2 


0.5 

0.5 

t3* 1 12x10"* 

2.0 

u 1 


1.1 

0.4 

2.0 


1.2 


0.5 

0.5 

14. t21x10 ‘ 

1.4 

1.2 


1.2 

0.5 

0.6 


0.3 


0.1 

0.1 

15. t31x1o’' 

1.7 

1.1 


0.1 

0.4 

0.3 


0.2 


0. 1 

0.1 

16, t4lx1o'* 

2.5 

1.3 


1.7 

0,5 

1.5 


1.1 


0.5 

0.5 

17. tsixio"* 

1.8 

1.2 


1.4 

0.5 

0.5 


0.3 


0.1 

0.1 

18. gjfxlO^ 

2.5 

1.3 


1.7 

0.5 

2.0 


1.2 


0.5 

0.5 

19. Tllxio""* 

62.0 

84.8 


35.3 

4.6 

22.9 


7.6 


1.2 

1.1 

20. Tl2z10”^ 

62.0 

84.8 


35.3 

4.6 

22.9 


7.6 


1.2 

1.1 

21, 121x10""* 

99.1 

41.2 


20.3 

4.1 

37.0 


12,0 


1.9 

1.6 

22. 131x10""* 

138.4 

28.3 


9.6: 

15.5 

108.1 


36.1 


5.5 

5.4 

23. 141x10""* 

49.5 

36.7 


16.8 

3.4 

20.0 


7.7 


1.3 

1.2 

24. 151x10""* 

179.8 

46.9 


30.8 , 

12.8 

111.3 


37.2 


5.6 

5.5 
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TiBIE 11 

Total Production Cost Minim zation ■ 
Unscheduled Tool-Eeplacenient Policj’’ 


TIPE-A, horse PO’®E=10~.o|AU¥diidADTc” TYPE-B, horse¥6^j¥r=io,o 

/ ■'rTTT-' At?, AATi;^TT^rv-■■rr'•^Tm ^ 


pTBEE OE COMPOHEIITS TO BE PBDUJ 

Iced = loooo 


P2TP1 
P:: I 

■>^0 Y 


= P,I 

M 4 a ^ X 


= 


IhUMBER OE GOT.CPOHENTS TO BE PRO DU- 

Iced = 200C0 


x11x10 

x 12 x 105 

x 13 x 105 

x 14 x 103 

x15x10^ 

x2lx10~ ' 

x22x10"^ 

x23x10”^ 

x24x10'^ 

x25x10“' 

. HPxlO” 

t 11 x 10“ 

1 12x10'' 

t21x10'' 

t31x10'' 

t41x10^ 

tsixio'' 

OSxlO^ 

T 11 x 10“ 

T 12x10"' 

T21X10"' 

131x10"^ 

T4lx10“_ 

T51x10" 

SLx10“^ 

PTx10”2 

PC 

TGC 

ITG 

TO 


2.2 

1.8 

7.5 

2.2 

75.3 

56.5 

61.3 
104.7 

87.5 

38.3 
2.9 
2.9 
2.8 
1.7 
3.2 

1.5 
3.2 

88.6 
88.6 

203.1 

294.1 

96.1 
133.0 

22,0 

35.2 
105.6 

30.6 

89.9 

261,5 


5.6 
2.1 
1.5 

6.7 

1.4 

78.1 

48.8 

54.4 
121.9 

92.1 
66,7 

2.8 
2,8 

3.4 

2.4 
3.0 

2.4 

3.4 

80.5 

80.5 

330.1 

496.2 

64.6 

161.4 

20.9 

37.2 

112.2 

36.2 

92.5 

278.3 


7.9 

2,4 

2.4 

6,0 

2.4 

51.5 

44.1 

57.1 

134.6 

201.9 

57.5 
3.2 

3.2 

3.3 
1.7 

3.1 
0.6 
3.3 

251.3 

251.3 

408.2 

289.6 

524.7 

9.2 

36.3 

108.9 

36.2 

57.4 

239.8 


8.3 

2.7 

2.5 
6.9 

2.8 
48.9 
43.0 
55,2 

137.5 

210.0 

61.7 

3.1 

3.1 

3.1 

1.6 
3.1 
0.5 
3.1 


36.1 
108.3 

36.2 
57.4 

238.0 


][=.5000' I 

12.8 

12.8 

12.8 

23.7 

12.8 
34.1 
' 34.1 

34.1 

34.1 

34.1 

28.4 

2.9 

2.9 
0.8 
0.4 
3.1 
0.7 
3.1 

464.8 

464.8 
751.3 

992.9 

540.2 

987.9 

142.1 

68.2 

204.6 

28.7 

16.8 

318.3 


Ip. = 


= R} = P3 I = P3 

n3^I=.9772^ |=.9987^ 


50.3 

50.3 

50.3 

121.9 

50.3 
27.9 
27.9 
27.9 

27.9 
27.9 

24.0 
0.9 
0.9 
0.2 
O.'l 
0.7 
0.2 
0.9 

312.3 

312.3 

504.7 

998.1 

296.9 
927.6 

90.1 
19.8 

59.4 
. 27.0 

0.0 

106.2 


131.8 
131., 8 
131.8 

178.7 

131.8 
25.2 
25.2 

25.2 
25.2 
25v2 
20,6 

0.4 

0.4 

0.1 

0.1 

0,6 

0.1 

0,6 

209.6 

209.6 

338.9 

989.9 

307.7 

981.7 

67.8 

13.2 
39.6 
27.0 

0.0 

79.2 


142.1 
142.1 
142.1 
181.0 
142.1 
24.1 
24.1 
24 -. 1 
24.1 

24.1 

20.2 
0.4 
0.4 
0.1 
o;i 
0.'6 
0.1 
0.6 

199.3 

199.3 

321.9 

97T.'5 

295.3 

985.3 
65,0 
13.2 
39.6 
27.0 

0.0 

79.8 
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Total Production- Cost Minimization - 
Scheduled Tool-E^lacement Policy (2 -a) 




TYPE-A H0E5E POWHitlO.OlilllOMJIIC lYIE-B 

— - .L i 1 -2 


1. xllxlO'^ 

2. x12x105 

3. x13x105 

4. x14x10^ 

5. xl5x10^ 

6. x2lx102^ 

7. x22x10""‘ 

8. x23x10~'^ 

9. x24x10“'> 

10. x25x10'“’' 

11. liPxlO"] 

12. tllxio' 

13. t 12x101 

14. t21x10'’ 

15. t31x10l 

16. t41x10^ 

17. t51x10^ 

18. 02x101 

19. Tl1x10”1 

20. T12x10“J 

21. T21x10~‘ 

22. T31x10"1 

23. T41x10~'^ 

24. T5lx10~1 

25. IS 

26. PTx10~2 

27. PC 

28. TGC1 

29. ITC1 

30. TCG2 

31. ITG2x10“1 

32. TGxlO-1 


77.0 

60.3 

62.4 

67.6 
62.4 

63.7 
3.0 

3.0 

3.1 

1.4 

7.4 
2.6 
7.4 

88,6 

86.6 

189.0 

236.6 

513.6 

H2.3 

34 

162,8 

488.4 
0.0 
0,0 

539.3 

118.8 

221.5 


f 1 = ^2 

ft 


4.2 

1.3 

1.5 

4.4 

1.3 
66.0 
64.8 
88.1 
74.6 

88.4 
333.3 

4.5 

4.5 
4.0 

1.4 

7.6 

2.6 

7.6 

169.7 

169.7 
192.9 

107.0 

412.5 

193.7 
6 

167.2 

501.6 
0.0 
0.0 

108.0 

20.5 

81.5 


Pi - Pgpi - P 2 |Pi = I 
xi772^^L9987^^ L 5000^ 


[.8413 


3.8 
1.1 
1.1 

4.4 

1.9 

71.2 

84.3 
110.9 

84.7 

55.4 

648.1 

4.5 

4.5 

3.8 

1.6 

6.7 

2.9 
6.7 

146.3 

146.3 
101.6 
68.1 

278.3 

618.3 
3 

148.4 

445.2 
0.0 
0.0 

54.0 

11.6 

61.5 


3.1 

4.2 

4.5 

6.0 

1.9 

4.2 

4.5 

6,0 

1.6 

4.2 

4.5 

6,0 

7.4 

6.0 

6.7 

7.3 

1.3 

4.2 

4.5 

6.0 

57.0 

40.3 

52.6 

61.1 

54.0 

40.3 

52.6 

61.1 

73.2 

40.3 

52.6 

61.1 

59.0 

40.3 

52.6 

61,1 

77.0 

40.3 

52.6 

61,1 

720.1 

635.6 

710.5 

806.7 

4.3 

7.4 

5.3 

3.4 

4.3 

7.4 

5.3 

3.4 

3.4 

2.1 

1.5 

1.0 

1.6 

1.1 

0.8 

0.5 

5.7 

10.3 

7.1 

5.6 

3.0 

1.8 

1.3 

0.8 

5.7 

10.3 

7.1 

5.6 

235.7 

635.9 

264.1 

133.3 

235.7 

635.9 

264.1 

133.3 

263.6 

981.2 

426,8 

215.4 

179.5 

983.3 

961.8 

629.4 

596.2 

896.9 

359.6 

211.6 

326.1 

997.3 

918.1 

647.7 

3 

3 

3 

2 

125.4 

226.6 

156.2 

123.2 

376.2 

679.8 

468.6 

369.6 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

54.0 

54.0 

54.0 

45.5 

11.6 

11.6 

11.6 

9.2 

54.6 

85.0 

63.9 

50.8 


HORSE PO\iiril=l 0 .d" 

k = i’2F^p7 

I = P 3 I = P, 
i.9772 ^f.9987 ^ 


6.5 

6.5 

6.5 

8.2 

6.5 
92.2 
92.2 

92.2 

92.2 

92.2 
905.2 

2.5 

2.5 
0,8 
0.4 

4.5 
0.7 
4.5 


2 

99.0 

297.0 

0.0 

0.0 

45.5 
9.2 

42.5 
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4.27, 4.29, 4.32, aad 4.34) with piobahilistic coefficients. 
the piobabilistic constraints for each cutting-tool are supposed 
to be satisfied at the same probability level, i.e., p^ = Pg = 

Pj* 5tour probability levels, .5000, .8413, .9772, and .9887, are 
selected for the present study. The probability level of .5000 
is equivalent to treat the constraint to be deterministic. 

The following conclusions, common to both objective 
functions and the the two tool- replacement procedures, can be 
drawn } 

1. The consideration of probabilistic tool-life has a signifi- 
cant effect on the machining parameters. The optimum 
cutting condition for the probability level ,5000 differs 
highly from the other levels under consideration, 

2. The cycle time decreases as the probability level increases, 
For example, in Table 9 the cycle time for the automatic 
Type-A goes down from .25 minutes to ,05 minutes as the 
probability level changes from ,5000 to .9987. 

3. Except a few cases, speed and feed rate increases with pro- 
bability level, lor exsmple, in Table 9 feed rate for first 
spindle for aitomatic type A goes up from .0079 ipr. to 
.0191 ipr. The same is true for speed. 

4. It can be observed by glancing throu^ iiie madiining times, 
t. ., in the Tables 9 - 12, that the variance in these time 
is more for automatic 1ype - B compared to the type - A. 
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It is becmse of availability of more freedom in. selection of 
machining parameters in Type - A. 

5, The cycle time for a given probability level is less for the 
automatic type - B as compared to the 1ype - A. The observa- 
tion remains unexplainable. 

6, The various tool— travel lengths for the operations to be per- 
formed on the job (Pig. 4) varies • ftom .125 inches to ,2.5 
inches. The drilling operation on the 4th spindle has got the 
maximum machining length, i.e., 2.5 inches. Therefore, it can 
be expected iiiat the 4th spindle would dictate the cycle time. 
The results compiled here are for 32 (4 probability levers x 2 
types of automatics x 4 tables for different objective functions) 
sets of values. Out of the total 32 sets, it is observed that 
the 4th station is having the maximum machining time 24 times. 

7, The results show that the madiining parameters are insensitive 
to -the edacity of the axtomatic. The Tables 9 and 10 are for 
aitomatic of 10 HP and 20 HP capacity respectively. It is 
observed that the values of machining parameters obtained for 
both the cases are the same. It is due to particular choice of 
job and cutting-force equation exponents. 

The conclusions dravn above show that a hi#ier pro- 
ductivity is possible in some cases with increased confidence 
in tool-life constraints. This also suggests that the current 
practices regarding selection of machining parameters based on 
deterministic tool-life values may yield a higher cost of pro- 
duction. However, an observation of opposite nature (a low 
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productivity with higher confidence level) is also possib;j_^ 
v/ith different set of data. 

For the sake of clarity, the results of the two 
objective functions are discussed separately. The 
of the results are given in the following sections, 

7.2 Cycle Time Minimization 

The Tables 9 end 10 depict the results for this 
objective function. The machining time is inversely 
tional to the speed and feed rate, aud therefore the Eiinijj^j^^^__ 
tion of qycle time is equivalent to maximizing speed end 
rate* The Table 9 shows that for high probability level 
(*9987), the speed for drilling operation on 4th ^indle 
very near to its maximum possible value (2400 rpm). ^ ‘^aage 
in ^eed from 1269 rpm. to 2399 ipi^* occurs for the drillsj^g 
operation as the probability level goes up from .5000 to ,993^^ 
Both types of the automatics reveal similar reailts for 
rate and speed, 

7.3 Total Production Cost lihniini nation 

This objective function is studied under two tyjj^g 
tool-replacement procedures. The Tables 11 and 12 depict 
results for unscheduled and scheduled tool-replacement 

cedure (2 -a) respectively. 

The magnitude of the various costs obtained are 
subject to time and cost paraasi®^® taken in this study. 
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However, it can be obsciwcd that the machine and labour cost 
is most important for the assumed input data. The second 
level of importance goes to the non— productive costs and 
obviously the least impoitant is the tool— cost, jin ext®!*^ 
sive sensitivity analysis on time and cost parameters is 
needed to draw any generalized conclusions rcgai’ding the re- 
lative importaice of these three costs. 

In the follovhng sections the two tool-replacement 
procedures are discussed. 

7»5»1 Unscheduled Tool-replacement Procedure 

It is observed in Table 11 that a very low cycle 
time is possible in aitomatic type-B compared to type-A. 

For example, for the various probability levels under consi- 
derations, the cycle times in minutes for ts^pe-A are ,32, 

• 34, .33> and ,31 vhilo for type-B they are .31, .09, .06, 
and .06, The author is unable to explain this observation. 

7,3,2 Scheduled Tool«!*replacement Procedure (2-A) 

The results for this procedure are presented in 
Table 12. A factor, 'C. is defined in section 4.2.2 of 
Ch^ter 4. It is written as 




vhere 


th .th . „ 

mean tool-life of i tool on 0 spindle, 
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^±0 ~ standard deviation of T.. , 

1 J 

SI - tool-replacement interval, and 

A full factorial experimental design was conducted 
for the factor q and the various probability levels under 
consideration. All integer values of q (-2 to 2) and the 
four probability levels (.5000, .8713, .9772, and .9987) were 
considered for the study. It was found that for a given value 
of probability level, the variation in q does not effect the 
optimal cutting conditions. This observation can be attribu- 
ted to the selected mean values and the standard deviation of 
the tool life exponents. 

The scheduled tool- replacement procedure is diarac- 
terizod by the duration of the period after ishidi all the tools 
are replaced. It is noted that the number of schedules required 
to produce 21000 (20000 fsultj^'w components + 5 % attrition 
allowance) components decreases from 34 to 6 as probability level 
is changed from .5000 to ,8413. This is also to be- noted iiiafe 
as probabili^ level increases further, the drop in number of 
schedules is not that steep, 

7.4 Comparison of the Two Tool-replacement Procedures 

The choice among the two tool-replacement procedures 
under study can not be decided in general, since the relative 
merit depend upon Uae various time and cost parameters, and the 
mean and the variances of the exponents in tool— life equation. 
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ho.vcvcr, "tile results obtained in liais study reveal that the 
totfl. cctt ui production is less in case of unscheduled 
t>^ol-rcplncuieu r .1 ccnpcjred to the sciieduled tool-replacement 
procedure. This be because of hi^h tool life vd.ues 
resulting in less tool failures aid also a high value (l25 
ninutus) chosen for the “tine required to change all ihe tools 
siruilt jnecusly for scheduled tool -replacement increases total cost 

i:«5 A 'tV: rd about Solution Procedure 


The aalution procedure does not guarantee global 
optinum. Therefore, the results obtained here may not correS' 
pond tc the global optimum. The optimum conditions can be 
iiaproved by starting ihe algorithm with different initial 
solutiv^n and thus selecting ihe minimum of all minimums. 

This improvement also cm not ensure to reach global optima 
but it vd.ll yield a better solution aid also the behaviour 
of ihe objective function will be revealed. 

The efficiesney of ihe algorithm can also be impro- 
ved by a proper selection of initial step size in Golden 
Seardi Section. The author has experienced that ihe initial 
value of step size affects greatly ihe total computer time 
required to reach ihe optimum solution. 

f,6 General CJonclusionB 

In this study ihe determination of optimum cutting 
conditions is integrated with tool-replacement procedures 
and ihey ceii be selected analytically* Puriheimore, it 



shows thet the probabilistic consideration of tool-life 
effects the OAtthag conditions, cycle time and, the total 
cost ui production significently. It can be concluded that 
the decisions based on detenainistic values of tool— life 
may yield a low productivity in scaae cases as observed in 
the present voik. 

As pointed out in Ch^ter 71, the data used for 
the analysis are primarily representative in nature. The 
involved arbitrari'ness in data prohibits one to make any 
generaliz at ionfi. 
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SCOPB OF EUTURE WORK 

8.1 Optiiaization Procedure 

A solution methodology can be used by an industry if it is 
economical and easy to adopt it. iny methodology suggested can turn 
out to be economicsl in two ways? l) the solution obtained by employ 
mg Hie suggested meHiod brings ^preciable reduction in cost of pro- 
duction and 2) the expenditures incurred in using the meiiiodology 
should be less ihan Hie saving in cost of production. The reduction 
in cost of production can be achieved by employing the algorithm 
suggested in the present work. Howeverj Ihe efficiency of the packgg 
be increcsod either by making iiiprovements in the algorithm or by 
using a different method of solution. The interior penalty function 
approach, thoi^h recommended by many authors, has been found sluggidi 
and time consuming in the present work. The author recommends to 
try anoiher method, such as Rosen's Gradient Projection method (40), 
to solve the problem. The relative merit of the two methods can be 
evaluated and the better of the two ^ould be selected. 

8.2 Optimum Grouping of Operations 

As mentioned in section 2,2 of Clhapter II, iiie area of 
optimum grouping of Hie ope rat ions in a multi-tool ease remains 
relatively unexplored. The problem also involves the redefining 
of operations and tool interference phencmena. The author feels 
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iaiat the graphical technique, referred as iSaalysis of Interconnec- 
ted Decisicn Areas (26) in the literature can be used to solve 
this problem. 

8.3 Application of Group Technology 

Aiother field open for fhrther research is the ^plicar- 
tion of the concepts of group technology in multi-spindle 
automatics. A lot of time is wasted in change-over of the pro- 
duction from one job to the other. The group technology concepts 
cm be used to group the jobs in such a way so as minimum altera- 
tions are needed for the change-overs. Once such group is found 
out, the technique of Branch and Bound can be put to use to 
determine an optimum seejueiiCe-of jobs to be manufactured. , 

8.4 Tool-life experiments for Ibim Tools 

The literature available shows that the form tools have 
received very little attentation. 'in experimental research can 
be conducted vhether the tool values for form tools differ sig- 
nificmtly from single point tool or not. 

8.5 T^fect of Horsepower of the Machine Tool 

In Ch^ter VII it was mentiesned, on the basis'of data 
taken for this study, -that the capacity of machine tool in 
terms of its horsepower, does not affect the results. This ob- 
servation does not seem to be convincing. Therefore, it may be 
worthwhile idea to investigate the affect of machine tool espar 
city using various sets of data. 
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iPPMDIX A 


DEEIYillOI OF METHOD OF CONYEHTIFG CHiUCE CONSTRAIITS IITO 
DETEIMIMISTIG C0N3TEAIMES 


Define random variables y^'s, as follows : 
n 

° 2 ’'j ■*’i ••• ’ “ 

3 * 


(A-1) 


viiere 


a. .'s are the technological coefficients, and b. Yg are the 

X3 X 


constraint coefficients. 


Since the sum of normal random variables is also 
variable, y^ is a normal random variable. The mean value and the 
standard deviation of y^ are given by 


y. 


n 

3=1 


a. . X. b. 
X3 3 X 


(A-2) 


cr 


^i t 


/ “ 

( 

''' 37k=1 


P o~ <3” X . X, 

' ^iJ. ^ik 3 ^ 


n ■ 2 N -1/2 

- !>’ p K O" rr-v X. +'y~', I 

a . b. a. . b. 3 b. / 

3=1 i3 X 13 X 1 ' 

Then the chance constraint (4.16) is written as 

/ ^ 2/ 

Prob (y.> 0 ) = / e""^ 2 dt p. 

/'■dO-"- 


(a-3) 


h(Pi) ,;v7f 


vhere 



Hence we have 
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end -thus 



bstituting equations (a- 2) aid (a- 3} into equation (a- 4), 
the equation (4.17), 



APPENDIX B 


CHOIOE OE INITIAIi VilDE OE ’r' 

To start the algorithm for unconstrained minimization 
a decision has to be taken about the initial values of 'r* and the 
factor c (the amount by #iich the value of ' r' is reduced at every 
step). The algorithm as such imposes no restriction vtiatsoever, on 
the value of and c. With large values of *r‘ however, the number 
of computations required to reach the constrained minimum increases 
tremendously, #iile.j small values of ‘r* distort the function 
bringing in hi^ eccentricity. The initial solution in the latter 
case has to be sufficiently close to the actual minimum. 

Eiacco and McCormick (l8.) have advanced two methods for the 

determination of an initial value of 'r^, were 

V 

V i.\) / j P ( )'l (B-i) 

vhere 

X is initial solution vector and 
o 

m _ 

P ^ ^0 ^ 1 / gj_ ( ) (B-2) 

The ecjration (B-2) chooses a value of r^ vhich makes the 
value of objective function aid penalty term equal. 

If r 0 , minimization of f (x) alone is carried out|, 
without considering the penalty tem. At every new point the value 
of r is checked by equation As soon as r > 0, the mini- 

mization of (x, r) is started with gradual decrement in r as 
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lias already been outlined. 


Wien r = 0, it injjlies that iiie unconstrained minima has 
been attained. 


Hoviever, if ihe starting point is close to the boundary 
the value of the penalty tem is excessive j thereby yielding a very 
small value of r. In such cases, the critical constraint is arbi- 
trarily assigned some minimum value. Thus 


f ( X ) / 


m 

Z! 

i=1 


mm 


TiTT 




^iiiia 


(B-2) 


Ihen the constraints are normalised any arbitrary value can 


be assigned to g^ , 



APPENDIX C 


CALCUIATION OE G-HADIENT VECTOR 

In order to calculate iiie direction used in Davldon - 
Eletctier - Powell's meHiod, it is necessary to calculate the gradient 
vector of the function f ( x ) during the course of each iteration. 

In -the absence of exact or explicit foimulae for the needed deilva- 

tives the gradient vector is ^proximated by the method of finite 
differences, 

xhe teciinique is simply to use the well known formulae of ■ 
differences to represent the derivatives. The simplest of Hiese is 
the approximation 

^ = ( f ( - f W / (0-1) 

■where 

X X^J . . . , X^ +II5 x^ , ,, .,x^) 

z = ( , x^, ...,x^) , and 

Ax^ = a small change in the variable x . 

3 2 

The following scheme is followed to calculate /\ x to 

3 

arit the problem studied in this -TOrk. The steps are ; 

1. Calculate A x. = vhere^is a sufficiently 

small daosen value, 

2. let xt = X. 4- Ax, , 

3 3 3 
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3. If = 0, modify A x to ^ is a 

very small number arbitrarily chosen, 

4. The constraint set is calculated and the xl is checked 

0 

for its feasibilily, if oonstraint is violated, 

gradient is dalculated by equation (c-l) and letting 

3 = 5 + 1 letum to Step 1, othervri.se go to Step 5'.' 

# 

5. If A X. is negative, modify A x. to - A x., fi md 

0 “ J ^ 

go to Step 4, otherwise Ax. = - Ax. and go to Step 4, 

3 3 

The value of and should be chosoi according to the 

function. In the present vroik ^ and S are taken to be as .1 md 

,01 respectively, Ihx (22) suggests that one way of ascertaining 

wiieliier the value of ^ is proper or not is to compute ' A ^ as in 

5 ^ j 

equation (C-l) fbr one value of Ax., recompute it for a smaller 

J 

value of it, md -tiien compare the results. If they agree sufficientlv 

tiien the ^proximation is probably sound, otherwise the values of 

•A X . need ad5ustment, Ibx has also suggested a ’ central formula’ 

J 

to reduce the 'truncation error’ involved in formula (c-l), but it 
increases tie mount of computation by 100 ;j. The 'central formula' is 


if f Cup - f 

viiere 


(c-2) 


= ( x^, Xg, ... , x^ -M x^ ) 

X. = ( x.|, Xg, ... , x^ -Ax^, ..., x^ ) 

However^ to avoid excessive computational effort, the 


first method is dhosen in liiLe preset study . 



appendix d 



COMPUTER PROGRAMME LISTING 





# ^ iz 5 *^ 5 ^ ^ 




LISTING OF T H £ S I S PROGRAM 

optimisation IN MULTI- SPINDLE AUTOMATICS 
-1 PROBABILISTIC APPROACH 

J}:* 5^ ❖=!' 4 V*^>!!*=ls5555!t^s^4!:4s* af sj: 5!= >!< :<t^5:^t ijtjjt 4^5|t 4;it :jc$3it4c :^f;5jt#* 5!!:i::{t *3i!**3S: *::Jt>[t 4!^>i'* = 

A’OMRd) -TYP; fJF STcOL 
A”’ATR(2) -OHN OF RAW MATERIAL 
AMaTRO) - YKLD STRcMGTh OF MATRIAL 
AMATR(A) - OlAM-rTRt OF THE RAW MATERIAL 
AMAF'US) - LENGTH OF BAR 
AMATRIE) - NO. UF SUCH BARS REQUIRED 
AMiATRIT) -RAW MATERIAL COST RATE 
CGMR(l) - COMPOliuNT RLCUIREMENT MATRIX 
CQflRd) - MU. UF CCMPO^eNTS REQUIRED 
COf'.HC:) - MATERIAL GRACE SPECIFICATION 
Cn.MRdn - CDMPf.iNtNT TOTAL LENGTH 
COHR(A) - material REQLIRcMENT 
CtTMR(5) - ATTRITION NO. 

C0MR(6) - OATt OF RECEIVING ORDER 
COMRIT) - OUc DmTS 


* 

3^„s§s, 3jt ^ 


COMR('j) - PRODUCTION STARTING TIME 
COMR(G) — GROUP NO. TC WHICH IT BELONGS 
COMR(IG) - FINISH TIME OF PRODUCTION 
COMR(il) -PRODUCTION PER UNIT 
COMR(l2) - COST OF PRODUCTION 

COMD(I,J,K) - DENOTES K DATA CHARACTERISTIC ON JTH OPERATION ON ITH STN. 

COMDCIfJ,!) - SURFACE PRODUCED 

COMDdfJiL) - TYPE OF CPMRATION 

COMD(I,J,3) - ACTUAL TRAVEL OF OPERATION 

C0MD(I,J,4) -DEPTH OF CUT 

COMUd » J»5) -TOOL USED 

C0MU(I,J,6) - FEED INDEX 

COMDd , J,7 ) - SPEED INDEX 

COMOdfJ»d) - MATERIAL REMOVING LENGTH 

C0M0d»J»9 ) - SECONDARY SURFACE OR PRIMARY SURFACE 

COMDd, J,iO)- WHETHER THE OPERATION REQUIRES DEPTH OF CUT 

COMDd, J, 11) - OPERATICN IS DONE ON AXIAL SLIDE OR CROSS SLIDE 

COMDd, J, a2) - SURFACE RMS VALUE 

COMDd,J,x3) - WHETHER SURFACE COMES UNDER SMOOTH CATEGORY OR ROUGH CATEGO 


C0MDd,J,A4) LENGTH TOLERANCE 
COMDd, J, 15) - OIAMETRE TOLERANCE 
COMDd ,J,i6) - HMAX 
COMDd, 0,17) - SPEED 


C0H0(I,J,13) - FEED 
COMDd, J, 19) - OLD OIAM. 

C0MD(I,J,2'J) - NEW OIAMETRE 

COMDd, 0,21)- SPINDLE REVOLUTION NEEDED 

COMOd ,0,22)-TIMe OF OFRATION 

CMCdjO) - 0 TH CHARACTERISTICS FOR ITH M/C 

CMCd,l) - NO.' -OF STATIONS . 

CMCCU,2) - NO. OF AXIAL SLIDES 
CMCd,3) -..TYPE. OF 'M/C LL ' ■' G. ' 



O c*>c^c^c~>c^c>c^cTf <nr <-> <r> <ni ir> c~><rs>c~>c~>c^c~5<r^c~ic~>c~>c^ c~> <r> c~*> <r> c~> <r> <n» <ni c~> <r5 c~» c~> <r^ <-> cr> <n( c~» 


c C:’X{It4) - 'jI*-F KL-'iT SFI jDLE RPM OR SAME 

I CHCCItE! '■iI '• ir'‘''.lM 'A/C RFM 

C CFL { i j6 } ~ ''A. si 'AJ‘'1 MAC F I '■i£ RPM 
cr-X(i,?) - Mir, MUM FtcX 
C'-’C{I,3) ’’AXiVij*' f.X-D 
CMC(i,S) - MJ.'iifUH DtPTH OF CUT 
CIX(1»1‘ ) lUf-i USPTh '.‘F CUT 

C'Xiijlr. UHRS.^ POWFR 
CFXCIjM) ‘!0. UF SPINDLE SPEEDS 
C:''X( I j !-;■) 'Q. JF FEEDS 

CrXd,!*-) RAPID TRaWcRSE SPEED 
CMC ( I j 1‘'. ) - PCW R TRANSMISSION EFFICIENCY 
CMCIMlo) SP.:l.I ['■■ if:D M/V PRODUCTION RATE 

CMCdlfdl - .dilMUM DIMENSION IN DIA. OF JOB WHIDH CAN BE HANDELED 
CfXIirlO MAX'. 1i)M DIM-KSiON OF JOB IN LENGTH WHICH CAN BE HACOELED 
CLXdd'-^) - AMAL MINIMUM FEED 
CtyCIid ) 'AAXIMUM AXIAL FEED 

TCMdiJjK) - i..'iiiJT-S KTH TOOL ON JTH OPERATION ON ITH STATION 
TCMd , 1 > - TYP,. OF TH^ TOOL 
TCMI (I » J» :) NO. OF TEETH 

TCMd»J»3) ACJY iJTHfR RELEVANT TOCL DIMENSION 
TCM(I,J»'d OI^METRE OF THE TOOL 
TCMd,J,5) RAKu ANGLE 
TCMd,Ji6) WEAR allowed 
TCM(I,J, 7} TOOL TIP RACIUS 
TCMd,J,b) MAXIMUM FORCE ALLOWABLE 
TCMd,J»9) ROUGHING TOCL OR FINISHING TOOL 
TCMd,J,l ) NO. OF RCGRINDS PERMISSIBLE 
TCMdjJfl’d TOOL EDGE COST 
TCMn,J,12) - HORCE POWER CONSTANT 
CuSTd) - DENaT.iS COST PARAMETRES 
COST (id DIRECT LABOR CCST/ 

COST! 2) OVER HEAD COST 


CUST(j) TOOL GRINDING COST 

C0ST(4) MACHINE SETTING COST 

COST{5) - INTrRML-DIATE ADJUSTMENT LABOR COST 

TlMECdtJ.K) DENOTES KTH TIME VARIABLE FOR JTH TOOL ON ITH STATION 

TIMECd.J/i ) SETTING TIME 

TIMECd,J,»2) GRINDING TIME 

TIMEC{I,J»3) IMTERMEDIATt ADJUSTMENT TIME 

TIMEC{,J,4} MACHINING TIME 

r Yf I T - CYfl F TTMF 

TOOLd,J,K) KTH PARAMETRE OF JTH TOOL 
AVERAGE FEED EXPONENTS 
ST. DEV. OF FEED EXPONENT 
average DEPTH OF CUT EXPONENT 
ST. DEV. OF DEPTH OF CUT EXPONENT 
AVERAGE SPEED EXPONENT 
- ST. DEVe OF SPEED EXPONENT 
TOOL (I, J» 7) - AVERAGE. VALUE OF CONSTANT IN TOOL LIFE EQUATION 
TOOL (I , J , « ) ST. , DE V. OF CONSTANT IN TCOLVLI FE EQUATION 
FORCEd.J.KJ - DENOTES KIH'PARAMETRE OF JTH TOOL ON ITH STATION 

FORCeIiIjIi) AVE.'-^yALUflf^^FE^^FONEWT. IN' force EQ. 


TOOLd f Jf i) 
TOOLd,J,2) 
TOOLdf J,3) 
TOOL ( I , J,4) 
TOOLd tJ»5) 
TOOLd ,J, 65 


.-i:; ' i' a.;xX-XX:.:,X'.' . v ' 



c 

c 

c 

c 

c 

c 

c 

c 


ra^C; { I , j,.. ) 31. n£v, f-ed exponent 
Fi;RC... (IjJjS) ^V.;. VALUEuFSP6EDEXPONFNT 
FORCr { I , J , S) ST. CcV. CF SPEED EXPONENT 
FJRC:{I,J,7) - --VE VALLE OF CONSTANT IN FORCE EQ. 

lilO-.-X?'’ - I'JtDc':* FOR MACFI.iE , WHETHR THE MACINE IS MULTI SINGLE TOOL 
IMDuXF - I‘'1CF.\ FDR OBJECTIVE FUNCTION 
IND'XC - I \Di,[ FUR UPrRA'flQIM TYPE 

I.iDTX - TtilS INUUEX IS FUR SOLUTION TO BE CARRIED TO THE NEXT PROBLEM 
Crj.'^MuN/Sr.RCH- /FUNVAL 

CCPI'^GN/RrPT/PU, lJ,TP,ITNCHtTEC,TGCH,GCH,SeTC,CIDLM,GRC,STEC,SGRC, 
ISSlTC* SCiDLM, STGCHfABC »NTEC , NGC , NGR , CCSTP , SGC6| NTNCH i AS YLT , FACTOR 
.If NS 

DIM^NSKirj PO( . 0),P{,iO) ,I YNCHdO) ,TEC( 10),TGCH(10) t GCH ( 10 ) , SETC { 10 ) 
},CIDLM{: )tGRC( ) fSTtCC..'. ) ,SGRC{lN),SSETC(ia),SCIDLM(in)fSTGCH(i:)) 
DfABCd )f';TEC(i ),NGC(1 ),NGR(1-) 

COyPCN/LIFE/TLIF ,TVARI ,FVARI,NTA 
CUMMUN/STS P/cTA 

COMHaN/TaR/FC...( . . ) ,Gi: (1 ),G22(l-: ) 

COW'MCN/TOTAL/NCONSTf f'iOej ,NTFN 
COMMON/PRil/PKfjd 
COMMON/'IIUI/R 
COMMON/CYC/CYCLT, AMAXO 
COMMON/INcULT/CON, SATS FY 
COMMON/MAlsE/NtNIN 

OIMLNSION NSPCJ{ .U) , AMA TR ( 10 ) , COMR( 10 ) ,COMD (8 , 6 , 25 ) i CMC ( 5, 2w ) t TCM ( 8 
1, 6f 25 ) fCOSKi ) jTIMFC ( 8,0,13) , TOOL? 8, 6, 15 ), FORCE ( 8, 6, 15 ) , PROB{ 10 ) , 
IT(1. flOfTCd ,NJ{2,3,5),C0N{l'5'J),HPi(iU,l()),HP2{lU,l0) 

DIMENSION X{5U) 

DIMENSION TLIFt (0) tTVARI (10) tFVARniC) 

DIMENSION Y( 5, ) 

COMMON/iNVAR/INXE, INYE 

common/di V/C l 

COMMON/MATV/AMATR,COMR ,COMD,CMC, TCM, cost, T I M tC, TOOL, force, NSPO 
COMMON /NVAR/NAMAT, NCOM R, NSTN , NOP T, NCOMD,NCMC , NTCM, NCOST ,NT IME ,NTP, 
1NFP,NMACIN,NSPNDL,NAXI L , NCROSS , N PROB , MN ,NSTI ON, NOP 
COMMON/INO/INDl:X,INDeXF,INDEXM 
COMMON/OUT/Y, lOY, I UX , N J, G1 , G2, FC ,HPl , HP2 , T, TC 
DATAC0STP/.5U/ 

DATA SGC6/50./ 

I NDEX0=4 
ETA=.25 


IND6XF=3 

READ i,INXE,INYE 
PRINT1,INXE,INY' 


PRINTl ,INXE, INYf. 

READ I ,NAMAT, NCOMR , NSTN, NOPT ,NCOMD, NCMC, NTCM, NCOST, NT IME ,NTPfNFP , 
IN MAC IN ,NSPNDL,NAXIL,NCROSS,NPROB 
PRINT I, NAM AT , NCOMR, NSTN, NOPT ,NCOMO,NCMC» NTCM, NCOST , NT I ME, NT P, 
1NFP,NMACIN,NSPNDL,NAXI L, NCROSS, NPR08 
READ 3, (AMATRCI ),I=1»NAMAT) 

PRINT 3,(AMATR(1),I=1»NAHAT) 

READ 4,{C0MR(I)tI*»l*NCCMft) 

COMRii J=C0MR{l)#2'l' “'TrYT V,. ;; ; 

PRINT A,(C0MR(U.*I=1*N50HRI^ . ^ 
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K 'AH I »(?'<'■)?)(*) f 1= . jr,(pR03) 

PHI-M 'Ij {PRUr.( I ) J I=i, , (ypi^Qg) 

4::A,'j ’3, i ) 5l = i rMSPNDL) 

PRIUT 5,('^SPC(I ) , 'gSPfNDL) 

5 FUR''’AT(., 15) 

na6 1 = . ».)Sp!’"jL 
\^i=jSPC{I) 

DIJ5J = .. ,,JN 

RLAD J,K),K=1,nC0MD) 

P R 1 1 IT 7 t ( 1-0 C f 1 , J , K ) , K =i , NCOMD ) 

I F (CllMC { I » J j4 ) . :q, 0. ) GO TO 6 

6 CUNTlOLr 
DOi)I = l ,'"1MACIN 

R::AD 9 , { C ( I f K ) , K = l , N CMC ) 

PRIOT 9,{CMC(I,K),K=*^CmC) 

8 COPiTINLE 
nuiOI=l,NSf’i'iuL 
MN = ilSPC(I 1 

DOi J=1,NN 

RtAD 11,{TCM( I, J,K ),K=1, =TCM) 
TCM(I,J,3)=TCM(I,JtB)/2. 

PRINT 1^, (TCMd ,J,K) ,K = i,NTCM) 

10 CONTI NLF 

RFAU 12,(CnST{niI=ltNC0ST) 

PRINT 12, (COSTd )»I=l,KCrjST) 

00131=1, NSPNDL 
NN=NSPC(I 1 
D013J=1,NN 

READ 1.4, ( TIMEC( I,J,K) , K=I,NTIME) 

PRINT 14, (TIH'.;C(I, J,K) ,K = l,NTIMe) 

13 CONTINUE 

D017I=1, NSPNDL 
NN=NSPC(I ) 

D017J=1,NN 

READ 14,(T00L{I , J,K) ,K=..,NTP) 

PRINT 14, <TOOL{ I,J,K) ,K=.,NTP) 

17 CONTINUE 

00191=1, NSPNDL 
NN=NSPC(I ) 

D019J=1,NN 

READ 14,(F0RCL(IfJ»K) ,K=i,NFP) 

PRINT 14, (FORCGd, J,K) ,K = 1,NFP) 

19 CONTINUE 

1 F0RMAT(16I5) 

2 F0RMAT(10F6.2) 

3 FORMAT (8F10. 2) 

4 FORMAT(2FlU.0,/l5F5.l) 

7 F0RMAT(8Fi0.3) 

9 F0RMAT(5F7«i,'itF10%'^»P5 »! ,/8F7 •lt2F8»4) 

11 FORHAT(10FG.35 , ^ 

12 ' FORMATI8Fi0.3 J ^ /.A'''' -O 
:14 . 'F0RMAT{8FlQi3) 



5- ' 
90 


5 


Go t:; ^5 ,& ,6*' ,6 ),indexo 

■STICK -I 

I'-.D \r-= , 

TL^)oT = ':. 

TTi 

■M?T = :-Si’3 C iSPlJlL) 

Di i" Ik= ,t <QPT 
\MP = I K 

CALL CCA5T-UX) 

P‘U iT 3‘”,A!rj,L 

CALL ‘ynip.Acx) 

PRINT ^ t I ’ D ; X L' 

PRINT ^ » C 0 S T { „ u ) , I N , N S T 1 0 N , T { N S P N D L , I N ) 

TLUST=TCaST+COST(lO) 

TTI.-L =TTI M:- + T {l-JSP'MnL, U) 

3i CONTINUE 

TlMcPJ=rTiyE/COMR( 1) 

COSTPJ=TCQST/CuyR{ 1 , ) 

STOP 

40 NSTI[jN = ,. 

GO TU 90 
50 N5TiaN=:. 
iNDf-xy = : 

TCCJST=0. 

TTIME=C. 

ihd:-xf=- 

D034I=1,NSPNDL 
NN=NSPC(I) 

D034J = i,Nr'l ■ 

CALL CCNSTR(X) 

PRINT 35tNINiN 
CALL yiNIMA(X ) 

PRINT 33,INDnXO 
PRINT 32rCnST{^,,), J,i,T(l»Jy 
TCOST=TCOST+COST( 10} 

TT1HE=TTIME+T(I , J) 

34 CONTINUE 

TIMEPJ=TTIME/COMR{ I ) 

COSTPJ=TCOST/COMR( 1) 

STOP 

60 D053MN*l»Hy<ACIN 
INDfcX=i 

NST10N=CMC{MN,l ) 

FACTOR=5. 

MCONST^O 
NOBJ=0 

NXFN=U, , , , , 

CMC(MN,6) =2400. CARE 

IF(CMC{MN,3).EQ.O.) CPC I MN, 11) =CHC (MN, 11 J /FLOAT! NSTION ) 

INOEXF = 2 ' " - ■''/ A', 

PRINT 33,IN0EX0 , A- :\.;:'‘:A' . ' 

CALL CGNSTRCX) , ' ' ' -a^ 



O O <-) 


6 


PRiniT 

D I. J • j T I • k = . » , • I "m 
39 Slll''i = 3UP + . ./Cu:'.( I.'O 

R = FUNVAL/5:.JM 
C... =-. 

CALL f'^rirntx) 

Du.; 

'ilM = ’'jSPC(I ) 

D U - / J " y J ' -I 

I F ( LliPC { 1 , .) t' 3 } . "Q, .) GO TO 37 
PPi-'iT 3C,CDST{„ )tJ»I»T(I,J) 

37 CUt'iTI.JL-; 

36 3 ON TI UL, 

PiU-'iT a?! iiCiJ.D'oTi!JQBJ,NTFi\ 

82 FO« ^AT ) 

I F( INDCXF.LO.C) GO TO 50 „ 

IF( iNUiXF. GO TO 5* ■ I 

PRINT 0 'fMCYCLT 
DOS 3II = ..tNTA 

503 PRINT 30 5 , TL I F ( 1 1) , TV ARl ( 1 1) , FVARl ( 1 1 ) , FC2 { 1 1 ) , Gil ( 1 1 ) , G22 ( M ) 

GU TO 5;. 

501 PRINT 5 '5,ASYLT,CYCLT 

504 D(:i6C6Ii = . ,HTa 

506 PRINT 503,TLIF{ II ) , TV ARI { 1 1 ) , FVARl ( I 1) , FCZII I ) , Gil ( 1 1 ) , G22 { 1 1 ) 
D05' 71 I=i ,NTA 

507 PRINT 5 5, TEC(II),TGCH(II),SETC(in,GRCtII),CIDLM(II),A8C{II), 
IITNCIK II) jNTECdl) ,NGR(I1) 

I F( INDEXF.cQ.3) GO TO 509 
GO TO 53 

502 GO TC 5-1 

509 0050811=*, NTA 

508 PRINT 505,STlC{II ),SGRC(II) ,STGCH(n),SSETC(II) ,SCIDLM(II ) 

PRINT 5".'P,NS 

50C0 FORMAT (IX, 511 .>) 

505 F0RMAT{IX,6F.12.2, 5110 ) 

53 CDNTINLE 

57 CONTiNLfc 

33 FORMAT (33X,40PERATI0N INDEX=*,I5) 

32 F0RMAT{1X,*TQTAL C0ST=*, F10.2,/IX,*TI ME OF OPERAT ION=«', 1 5, 2X, *0N*, I 
I5,2X,4STATI0N4, F10.2) 

35 FORMAT { 10X,*N0. OF C0NSTRAINTS=*,I5,/10X,*N0. OF VARIABLES=*,I5) 
STOP 
END 


SUBROUTINE INTAL(X,Y) 

THIS SUBROUTINE INITIALIZE THE VARIABLES 
COMMON/iMVAR/lHXe, INYE 



D1 ) 

Di SI!. -i Y{ir, ) 

X ( „ ) = .i: 6 

Xil 7 

X C / 

/■; {:)-=. 

Y( ?! = :': . 

Y( ..)=6''' . 

Y{.'. )=oC.’. 

Y(‘,-)=dY- . 

Y(5)=6': . 

Y{i)=^': . 

i") t j J. I = * I i ''■I l" 

1 / ( I 1=ALQS{X( I ) ) 

Du:- 1 = - »iny:, 

2 Y( 1 )=ALOG{ Y( i) ) 
R^;TUR^J 

c ("1 D 


SUBRQUTIi^t: QBJFR(TCQST) 

C 

C TIS SLBROUTIME EVALUARES THE TOTAL COST FOR TWO SCHEDULING POLICIES 

C 

COPMON/KEPT/PU,P,TP,ITKCH,TEC,TGCHtGCH,SeTC,CIDLM,GRC,STEC,SGRC» 
iSSETC»SClDLM,STGCH,ABC ,NTEC, NGC, NCR, C0STP,SGC6tNTNCH,ASYLT, FACTOR 
ItN.S 

CUMM0N/LIF::/TLIF,TVARI ,FVARI,NTA 

DIMENSION TLIFUff ) ,TVARI (IC ) » FVARl ( 10 ) , FXX ( 30 ) , POI 10) , P { 10 ) , ITNCHI iO) , 
I( iO) ,TEC(iO),TGCH{ 10 ) , GCH( 10 ) ,SETC ( 10 ) ,C IDLMI 10 ) , GRC ( 10 ) , STEC J 10 ) , 

ISGRC (1C ) ,SSETC( IC) , SCI CLM ( IG ) , ST6CH( 10 ) , ABC ( 10 ) , NTEC ( 1C ) »NGC ( 1C- ) , 
iNGRdC ) 

COMMON/INVAR/INXE, INYE 
COMMON/CYC/CYCLT,AMAXO 
COMMON/MANE/N,NIN 

COMMON/MATV/AMATRfCOMR,CQMD, CMC, TCMt COST, TIMEC, TOOL t FORCE, NSPO 
common/nvar/namat,ncomr,nstn,nopt,ncomd,ncmc,ntcm,ncost,ntime,ntp, 

1NFP,NMACIN,NSPNDL,NAXI L, NCR0SS,NPR08, MN, NSTION,NOP 
COMMON/IND/IND£X,INDEXF, INDEXM 
CDMMDN/OUT/Y, lOY, I OX ,N J,G1 ,G2, FC, HPl, HP2 ,T, TC 
COMMON/INEOLT/CON,SATSFY 

DIMENSION NSPOdO) ,AMATRdO) ,COMRdO) , COMO ( 8,6,25 ), CMC ( 5 , 20) , 
iTCM(8,6,25 ) ,COST (10 ) ,T IMEC (8 , 6» 10 ) , TOOL( 8, 6, 15) , FORCE( 8, 6, 15 ) , 
2PRO8dC),C0N(leO),Y(5id,HPld0,l0),HP2dO,10),X(50),T(i0,10),TC{iO,lO) 
3,i0) ,NJ(3,8,5) - ■ ■ 

ASYL=0. ' ' ^ 



2: 

1 


iw 


5 
4 

6 

20 


K '] = ..■ 

D’J. 1 = . ,^;ST1CN 
'4;;=''.SPC{i ) 
do; j = .„ 

iF(C'JNi3( I, J,':)) .cO.O. ) GJ TO 2 
K'i = KO+.l 

«GYL=;SYL+, ./TLIFIKN) 

AGYLT=1./ASYL 

TPR0T = 1.1*CYClT*CQMP.(} ) 

•'jTiaCli = TPR'IT/A3YLT + 3,. 

PRrjC = 7PRaT=!‘CGST (^) 

GO TC ij ,0 ),If\D£XF 
K0= 

DD*^! = „ tOSTIUK 
M''1=JSPC{I) 

D05J=»tN0 

I FCCOYCd , Jt9).cQ..: . ) GO TO 5 
KN=KN+i 

TPROT=T(I f J)*COMR( ) ) 

I Ti\CH(KN)=TPKOT/TLIF(KN)+i.-i. 

AINU=1. 

ABC{KN)=TPR0T/TLIF{KN) 

IF (A3C (KN) .LT.TCM{ I, J, 10) 3AIN0=0. 

MTECIKN ) = ABC {KN3/(TCMU, JfI0)+l. ) + l. 

TtC{KN)=FLOAT (NTECIKN) )*tCM( If J,ll) 

NGRIKN )=ADC( KN)-TCM{ I , J, 4.0)+i.-l. 

GRC (KN) =FL0AT (NGR( KN) ) *CCST(3)*TIMEC( I, J,2)*AIND 
SETC{KM = FLQAT(ITNCH{KN) )*TIMEC( I, J,3,)*C0ST(5) 

CIULK(Kf))=(CCSTP+C0ST(7) ) *{ T IMEC ( I , J, 3>/CYCLT ) *FLOAT ( ITNCHI KN ) 1 
TGCH{KN)=0. 

CONTINLO 

CONTINUE 

TCOST=fROC 

D06I=lfKN 

TCOST=TCOST+TEC ( I ) +TGCF ( I }+GRC( I J+SETCI I )+CIDLM{ I ) 

CONTINUE 

RETURN 

NS = i.i>i'COMR{l)/eXP(Y{l)) 

KN=u 


D07I=i,f'4STI0N 

NN=NSPOd) 

D08J=ItNN 

IF(COHC(I f Jf9).EQ.O.) GO TO 8 

KMsKN+l 

P0{KN)=1.-{TLIF{KN)-P(KN)*TVAR|{KN)/FACT0R)/(Y{1)#CYCLT) 
IF{PO(KN) .LT^O. 3 PO(KN}»e. 

IF(PC(KN3 .GT*l. 3 P0(KN3-1. 

ITNCHIKM) s=P0(KN 3,*FL0AT (NS)+l»-3.» , 

AIND=1. ^ ^ 

IF(FLOAT(ITNCH{KN3Kii.||5B*JF3^^^ 

NTEC ( KK) = { P0( KNl*Fi,OAT-;|^l,:5./*T€;>tirt4f lQi+1- 



J.il ) 


} A 1 N C 


*’4 OR ( KN ) = FL'jflT ! I TmCH { KM :) ) -TCH{ I » J ,10 ) + 

T..C{Kr\)=FLGAT( iT..C{KN) )«]CM{I,J,il) 

STACtKN ) = FLtjA7 ^^5)/(TC^{ I } + *, l^TCMC I • 

STGCH( KN) =f'. 

GRC (KK )=FL 0 # 1 T (!^';R{ KN) ) = 4 CCST ( 3 ) IMEC ( I, J, 2 )*A 1 NC 
SGRC(t<M={FLCaT( jS )-TC^(i, J, 5 , ) ) =«COST (3 ) *TIHtC <1 , J , 

Sc TC { K ^ ) = FLC T ( J T 'JCH { K N ) )*TIN';£C ( I , J , 3 )*COST { 5 ) 

SScTC ( KN) =FLC^T (fsS )*TI NEC( I, J ,3 )= 4 =CCST( 5 ) 

CIDLFCKN) =(Tiy,_C(I ,J, ^ )/tYCLT )*FLOAT( ITNCHCKN) )*(CQSTP + CG: 
SC IDLF { KN ) = { ClSTR + COST (7 ) ) =kfL 0 AT (NS ) *(TIMEC ( I , J ,3 ) /C YClT ) 
TCCH( Kh )= . . 

8 CuNTIKL'f 

7 CUNTiM.t 

SC(.JST = SGC&=KFLCAT{MS )*CCST(3) 

TtCST=C. 

0091 =,. ,KN 

TUCST = TUCST + T::C(I) +TGCF{ I) 4 GRC{i ) + S£TC{ n+CIDLMti ) 

9 CGNTIMF 
TSCST=SCOST 
DCl,iI = l,KN 

TSCST = TSCST+ST:iC(I )4STGCH( I )4SGRC ( I ) 4SSETC ( 1) 4SC IDL w { l ) 

U CGNTIME 

TCCST=TSCST4TLCST4PRCC 
30 RfcTURN 

END 


SUBRCITINE CLTPUKX) 

C 

C THIS SUBROUTINE PRINTS VALUES OF DIFFERENT VARIABLES 
C ■ : 

COf'NON/INVAR/INXE, INYE 
COPRGN/MANE/N,NIN 

COMMON / MA T V/ 4 FAT R, COMR ,CC!YC, CMC, TCM, COST, I IMcC, TOOL, FORCE, NS PC 
COMMDN/NVAR/NAMAT,NCOFP,NSTN,NOPT,NCOFD,NCMC,NTCM,NCCST,NTiME , 
lNTP,NFF,NMACIN,NSPNCL,NAXIL,NCRCSS,NPROB,MN,NSTION,NCP 
COMMCN/IND/INCfcX,INDcXF, INDEXF 
COMMON/OUT/Y,ICY,IOX,NJ,G1,G2,FC,HP1, FP2,T,TC 
DIMENSION NSPC(ifi) , AF^TR ( IQ ) ,COMR( 10 ) ,COMC ( 8 , 6,25 } , CMC ( 5 , 2u ) , 
1TCM<8,6,25) ,CCST(lN),TIMcCl8,6,lC),TCCL{8,6,15),FGRCE(S,6,l5 ), 
2PR0B{IC) ,C0N(1CC), Y{50 ),HPX(iO,lC) ,FP2C10 tIO),X( 50) ,T(lC,iC),TC( 10, 
3,10) ,NJ(3, 8,5 ) 

DIMENSION XX(5J) ,YZ(5a ) 

GO TC (10,20) iINDEXF 
20 DC1I=1,INXE 

1 XX{I)=EXP{X{I )) 

IDX=IN>E+1 

D02I=ICX,N 

I Y=I-INXE 

2 XX{I)=EXP{Y(IY) ) 

PRINT 3, { XX(n',I=l ,N) 


Ca 



F QRf^ AT { X , TPUT* , 5F 1 5. 4) 

RuTURa 

END 


CARE 


3 


SU.^RCLT1 'i,- C’-JiSTR(X) 

C ijc Jj! 5l<4e5t« 4>!t *!jei!c***si5 a}! 4:$ 4: 3jts5c;j.^aj ;{s :{(:{:){! 4 :$ :}!;}! :{e^sj!:)t:{c 

c THIS SLBROUTlNfc IS TO GENERATE CONSTRAINTS FOR A MULTI-SPINDLE AUTOMATOC 

c Machine having independent drives for each spindle 

dimension NSPGC , AMATRdO) tCOMRdC) ,COMD{8t 6t25 } »CMC(5 t2U) » TCM( 8 
1,6,?.5) ,C'JST(i ) fTIMECt 8,6tI0) »T00L(8»6,15),F0RCE{8»6tl5) fPROBtlO) 

DIMENSION CON d *0) ,Y( 50 tHPK 10,10 ) ,HP2 (10, 10), X{ 50) 

DIMEKSIOiJ T( iU, 10) ,TC{ 10,10) 

DIMcNSION NJ(3,8,5 ) 

COMMON/TOTAL/rJCONST,NOEJ,NTFN 
CDMMQN/SERCHZ/FUNVAL 
COMWQ J/CYC/CYCLT,AMAXO 
COMMON/PRQ/PRCB 

COMMON /MA TV/A MATR, COMR ,COMD, CMC, TCM, COST, TIMEC, TOOL, FORCE, NSPO 

ccimmon/nvar/namat, ncqmp,nstn,nopt,ncomd,ncmc,ntcm,ncqst,ntime,ntp, 

INFP ,NMACI N,NSPNDL,NAXI L,NeROSS,NPROB, «N,NSTION,NOP 
COMMON/IND/INDEX,INDEXF,INDEXM 
COMMCN/OUT/Y,IOY, I0X,NJ,G1,G2,FC,HP1,HP2,T,TC 
C0MM0N/SUB4/SIGMA 
COMMON/MANe/N,NIN 
COMMON/INEULT/CON,SATSFY 
CDMMON/INVAR/INXe, INYE 
DIMENSION PRCBHPI 5) 

COMMON/LlFE/TLIF,TVARI,FVARI,NTA 
COMMON/TOR/ FC2 ( lu ),G U (lO) ,G22( 10) 

DIMENSION TLIF{ 10) ,TVAPlIlf3) fFVARKlO) 

NTA-0 

IF( INOEX.EQ.O.ANU. INDEXM.EQ.2) GO TO 25 
IF(INDEX.EQ.l.AND.INDEXM.EQ.l) GO TO 31 
IF( INDEX. cQ.O. AND. INDEXM.EQ.i) GO TO 400 
j CALL INTAL{X,Y) 

!C INITIALISATION OF X ANC Y 

INDeX=C 
GO TO 27 

31 CALL INITtX) 

INOrX=0 
GO TO 27 

25 D028I=l,I0Y 

IDX=I0X+I 
Y(I)=X{IDX) 

28 CONTINUE 

27 iFdNDEXM.EQ.i) GO TO 400 
IOAXIL=0 
NCONST = NCONST+l 

NCo,N=c= ^ v. i,'.; 

NSpNDi=X 



12 I=uS)^i\CL 

=!\SPC(I ) 

PkubHP (.. ) =P. 

SHD = v', 

prOhmp , 

PP'jbHO {i;-) = . 

PtUJbHP { '.) --0. 

r.iij,; J = , j NN 

irccclwcd, J,9).,.G. .) GO TO 2 
IF(C0-''^C(I ) GO TO 3 

IF(CO--'C(I , J,6).GT.:. ) GO TO 4 

‘<JT/>=f>iTA+: 

irix=i + ':-iuAxiL 

PHlPl=FRUu{l) 

P(1I PI = PP'Jrt( 2 ) 

CHP^TC^d , 

N|Jdtl,J)=TNX 
IUX = IfO< 

NC(i.n1 = NCOU + ,. 

c MAXIMUM FiPD CONSTRAINT 

CC!N{NCCN)=ALGG{CMC (MNf 6) }-X( INX) 
NCON=KCON+" 

C MINIMUM FbLD CONSTRAINT 

:ON(NCCN) =-ALQG(CNC(MN,7 ) )+XdNX) 

GO TC 5 

4 INX=10X 
NJUd » J)=njx 
NTA=NTA+1 

5 IFCCOMCd ,J,7).GT.l.) GO TO 6 

INY=I+l 

N J(2,I » J)=INY 

IOY=INV 

NCON=NCaN+i 

C MAXIMUf' SPFED CONSTRAINT 

C(JN{NCCN) =ALOG{CMC (MN, 6) )-Y{ INY) 

NCON=NCON+l 

C MINIMUM SPEED CONSTRAINY 

CON(NCON)=-ALOG(CMC( PN,5) S+YdNY) 
GO TO 8 

6 INY=IOY 
NJ(2,I ,J)=INY 
GO TO 8 

3 INX=1 

NTA=NTA+1 

iaAXlL=IOAXlL+l 

NCON=NCON+l 

C MAXIMUM FEED CONSTRAINT 

CON (NCON) =ALOG(CMC (MNj 20 3 )-XCl) 

NCON=NCON+l 

C MINIMUM FEED CONSTRAINT 

CON(NCCN) «=-AL0G(C«C(HN,l9) )+X(l) 
IF(CONDCI t Jf7}*6T.l- ) GO TO 8 
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I')Y = I + 1 
N J ( .■'■ » I , J ) = I ■'iY 
i:jY = I'\Y 

"ICO = , 

M.'XU'Uh' 5 P'’..C CnjjSTRAIM 
cuia:icc:.i) ^ 1 LCG{CMC{MN, 6 ) )-Y(INY) 

'■■JCi,lu' = i\CLif'i+ . ■ 

MI.'jlYUY SR.£D CONSTRAINY 

C;jl. {MCC"'l) --=~ALGG{CHC{MN ,5 ) )+Y( INY) 

MCOM = NCiJM+ , 

G„ = TCCLII f J,& )** 2 + (TOOK I , J ,6 J ^^'Y { I NY ) )** 2 +( TOOL ( I ,J, 23 *XnNX) )** 2 +lTOOL( I, 

1( TGlIK 

Gi, = rrcCL{ 1 , J, !)*AL 0 G{ 1 C. ) )*=«= 2 -T 00 L( I, J, 8 )** 2 +G 1 
TOOL LIFL constraint 

CUNCNCCN) =CLOG( 1 I’l'TOCLd, J, 7 ) + M.+TCOL( )*XUNX) + (l.+TOOL(I 

i, J, 5 ) )«Y{ IilY)+TnOL(I ,J , 3 )*C 0 MD{I ,J, 4 )-Y( 1 )+PHIP 1 *(SQRT{G 1 } ) 

TLJ F:: = COrUMCQN)-PHlPl=«^(SCRT(Gl) 3 
TL1FC = TLIF!:+Y(1I 

CON( KCCiU=CUN(NCON 3 -ALCG(GOMD{I, J ,8 3 3 
TLIF{NTA 3 =FXP{TLIFc 3 

G 2 = F 0 RC:.:( I, J, B) ^* 2 +{FORCe( I, J, 63 *Y( INY) 3 ** 2 +( FORCE! I, J, 23 *X(INX 3 3**2 
1 *. + (FORC£ ( I , Jt 43 *COMO( I, J ,43 3 Y ’!'2 
G 2 =G 2 -FGRCE { I f J ,8 3 *Y 2 +{F 0 RCE{ I, J ,8 3 *ALOG(iO. 3 3**2 
\ICON = NCON + 3 

MAXIMUf^ FORCE CONSTRAINT 

GUN{NCCN3 =ALOG(TCM{I, J ,8 3 3-AL0G(10.3*F0RCE( I,Jt7 3-FORC£( I,Jj5 3*Y{IN 
lNY5-FDRCc { I , J ,i 3 *X ( I NX )-FORCE ( I , J, 3 3*COMD ( I , J,4 3+PHIP2* ( SORT ( 62 3 3 
FC=( 10.**F(iRCu( I,J, 7 3 3 *(tXP(Y(INY3 3 ** FORCE ( I , J, 5 3 3*( EXP (X { INX 3 3** 

IFORCE! I,J, 1 3 3 *(EXP(COMC{ 1 , J, 43 3 **FORCE(I, J ,33 3 

HPl { I ,J 3 =C! 4 P*FC*EXP(Y tlNY 3 3 *{CCMD( I ,J ,19 3 +COMD( I, J ,20 5 3 / 2 . 

TVARI {NTA}=TLIF(NTA)*SGRT(G 1 ) 

FVARI {NTA)=FC*SQRT (G 2 ) 

G 2 Z{NT 4 )=G 2 
GlI{NTA 3 =Gi 
FC 2 !NTA)=FC 
NCON = NCrjN + 3 

SURFACE ACCURACY CONSTRAINT 

I F(COHC( I ,J , 13 ) .EQ.O . 3 CON( NOON ) =-X( I )+ALOG( 2 .*COMO ( I , J , 16 3 
1 / ISIN{ 2 .*TCM( I, J, 5 3 3 3 3 

iF(CDND{ I, J, 13 ).EQ.l. ) CON (NCON )= 2 .*X{ 1 3 +ALOG(COMD( I , J, 163 * 8 . 

1 *TCM( I,J, M) ) 

PR 08 HP (I) =PR 0 BHP {1 )+ALC 6 ( 10 . 3 *HPl( I , J ) 

PROBHPd 3 = ( PROBHP( i)*FCRCEl I, J» 8 3 3**2 


SHP=SHP+PR08HP{ i)*PHIP2 
K=NJ{2,I,J) 

PR0BHP(2) =PROBHP(2 5+Y{K)*HPl( If U ) 
PR08HP(2) -( PROBHPC 2 3*FCRCE{ I, Ji 6) 3**2 
SHP=SHP+PR0BHP{23*PHIF2 
K=NJ(1,I ,J) . ' 

PROBHPO 3 ^PRpBHMS l+XI K| WPin» J 3 
PROBHP C3) =( 2)J**Z ^ 

SHP*SHP+PROBHPt31*Ptt" 

^ P RO BHP ( 4 1 * PRPPpfti 4l:t« 



I f J 3 
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P=’0!3HP(s-) = (PR.lBHP{-f)*FCRCE{ I, J, 4) )’«'=«'£ 
SHP=S^iP+PRCBHP{ 4)-'PHIP2 
SHP = SuRT( SMP ) =;-3(^RT (PHI F2 } 

HP: ( I s J) = .iHP 
2 CDNTit'jf;- 

N = ;jSPU{I ) 

D0_ J = l,li 

iF(CCHC(I , J,9) .3Q. .) GO TO IQ 

NCO'\ = NCQN + ; 

COi'HfNC'J 'j) =CKC ( KN, il ) -HPK I , J )+HP2( I »J) 
C HORCi; PUWv;.R CCiJSJRAINT 

iJ COMTIMJc 

^SPNDL = NSP'lDL + _ 

I F(NSPNDL.GT.fiSTION) GC TO 11 
GfJ TU 12 
11 D026I=l,inY 
IUX=IGX+I 
X ( IUX}=Y{ I ) 

26 COMTINt.:: 

FUNVAL=FUNXON( X) 

DD45I=1,MSTIDN 

MN=NSPC(1) 

D046J=1,MN 

IF{COf^D( I , J,q) .cQ.' . ) GO TO 46 
MCOH=NCf.lN+:. 

CO^vi{NCUN)=CYCLT-T ( I, J )+.*>l 
46 CONTINL2 
45 CONTINUE 

NCDN=NCOfl + „ 

C NO. OF COMPONcNTS CONSTRAINT 
CON{NCOM)=ALOG(COMR{l ))-Y(l) 

NIN=NCCN 

N=IOY+IDX 

RETURN 

400 I=NSTION 
J=NOP 

PHIPi=PROR(l) 

PHIP2=PR08{2) 

NC0N=O 


NN=5 

I F(COMC(I , J,9).cQ.O. ) RETURN 
CHP=TCR(I , J,12) J 

NCON=NCON+i 


C MAXIMUM FEED CONSTRAINT 

GON(NCCN)=ALOG(CMC (MN,8) )-X(l) 

NC0N=NC0N+1 

C MINIMUM FEED CONSTRAINT 

C0N{NCCN)=-AL0G(CMC(MN,7) 5+X{l) 
NCON=NCON.+i 

C MAXIMUM SPEED CONSTRAINT 

C0N{NC.CN)«AL0G( CMC {MNf 6J I-XUJ 
3' '' . NCON=NCON->-l . ■ y 

c MINIMUM SPEED CONSTRAINT ; 
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ClI‘‘ 1{KCG j)=-'\LDG(Ci«C( )+X(2)' 

Al Ct JN = W C lJ''J + . . 

c U.-.PTH GF CUT CCNSTRAINT 

CON{-GCC -!)=AL0G(CMC(MM,I. ) )-X(3) 

'MCUi‘j=NC'j. j + ' 

C 3-;'’TH OF CUT CCNSTRAINT 

CUh{ I iCCi'I) =-AL0G(CMC(MN ,9 ) )+X( 3 ) 

NCn,'j = NCnN+ , 

G ;=T(:CL( I t J ,3 )*-'2+{TQ0L( I , J,6)«x (2) )** 2 +{T 0 QL{ I, J,2)=!‘X(l) )**2 

l+{TC!jl_{i}Jt4)^X!3))4#2 

G:>FORC.;.{I»J»8)+{FaRC£{l,J,6)=«'X(2) )**2+{FORCE{ I,J,2)4X(1) 

12+{FORC -{ I , )4’i'2 

C TOOL LIFE COnSTRAI nT 

COOdlCCN) =-A lCG(C0MD( I ,J,3) )+AL06(i0. )*TGOL( 1 1 J »7 ) + ( i.+T0QL( I, Jt5 ) 

1)4X(2) + (1. + T0CL{1, J,i ) )*X(i)+T00L{I,J,3)’!=X(3)-AL0G(X(4) ) +PHI Pl=i‘ { SQRT ( Gi ) ) 
?.PT(G'.. ) ) 

CQMiNCC ■!) =CfJN(NCON)+ALCG(COMD{ It J»3 ) )-ALOG(COMD{ I »Jt3) ) 

NC(jN = NCOM + :, 

C MAXINUP FORCfc CONSTRAINT 

CON(NCCN) =ALOG{TCM( ItJ tS ) )-ALOG( IP . ) *FORCE ( 1 1 J 1 7 )- FORCE ( I , J , 1 )*X { 2 ) 
D-FGRCE (I ,J,5 )4X(1 )-FORCcCl f Jt3)*X(3)+PHIP2*{SQRT{ G2) ) 

FC={10.*4FORCr:{ I ,J ,7) ) =»{cXP{X( 2) )**FORCe(I , J, 5) )4(eXP(X(l) )*4 
IFORCEdtJti ) )4(t;XP{X{3))’!'*F0RCE{I,J,3)) 

HPK I, J) = AL0G{CHP)+AL0G(FC)+X(2)+AL0GtC0HD(It4t7)+C0MDl It Jt8)/2. ) 

HP2(It J) = PHIP14(liPl(I t J) )*{ (FORCE! I, Ji8) )=i<*2) + (X{ DEFORCE (1,4,2) )**2 
l*2+(X(2)*FORC£( 1,4,4) )’!*2+(FQRCE(I,4,6)=i'X(3) )*42 
NCON = NCON + 4. 

C SURFACE ACCURACY CONSTRAINT 

IF(COP>C{I ,4,13) .EO. ■ . ) CON(NCON)=-X{i>+ALOG(2.*COMD(I,4,l6) 

1/ (SIW{2.*TCM(1,4,5 ) ) ) ) 

IF(COPD( I,4,13).eQ.i. ) CON(NCON)=2.#X(l) + ALOG(COMD{ I,4,16)*8. 

1#TCM(I ,4,M) ) 

NCON=NCON+i 

C NO. OF COMPUNcNTS CONSTRAINT 

CGN(NCCN)=ALOG(COMR(i) )-X(4) 

NCON=NCON+l 

C TOOL VIBRATION CONSTRAINT 

CON(NCCN)=ALaG{TCM( 1,4,3 ))-X(3)-X(l) 

NCON=NCON+l 

C HORCE POWER CONSTRAINT 

CON ( NCQN ) =CP!C( MN, II ) *CMC(MN, 15 )-HPl( I, 4 )-HP2( 1,4) 

N4{l,I,4)=i 

N4(2,I,4)=2 

N4(3,I,4)=3 

NCON=NCON+i 

CON(NCON)=X(1) 

MCON=NCON+l 

CON(NCCM)=X(2) 

: NCON=NCON+l 

CON{NCON)^X{3) 

S' , NCON=NCON+l , 

CON(NCCN)pX{'<») 

'N= 4.';v,. „ r 





? 5 

i.ii H=.‘'.CC'''; 


S'J'i'-ll'L'TM - CQ ,STR( X) 

Q '4- '4‘ 4* 4< - 1 « ❖ 4« ^ ^ -v ^ ^ ^ 5}e ijt jjt jjc ^c j§c .5^ 5^ 4 « .^c 4c:4t4c5jc4«-45 4c = 1 ^ ^ ^ 

Q4:4<4f^« 5:fi4«5^^4 4 4< 4« 4< 4« ❖ 4' 4^ 4^ 4- 4' 4c 4s 4 4* ❖ 4c 3^ ❖ 4c 4 e 4c 4c 4< 4 ? 4? 4c 4c 4c ^c .^c 4 c 4c 4 c 4 : 4 c 4 ; 4 : 4 c 4 c 4c' 4c 4c 4= 4s 4c 4c 4« 5^ 4s# 4c 4c 4s 4c 4c 4= 4c 4s, 

Diy^.NSrjri :<{2 ),Y(Z ),NSpOav),AMATR{ l0),C0MR(2<)»C0MD{5,2,25)t 
jC'|C< 2,2 ) ,TCR{5,2,?.5),C0ST(i..),TIMeC(5»2»i3)»T00L(5,2,l5),F0RCE(5,2,i.5J 
j2,_S) ,Pk:pU„G),T{] 0,10 ),TC( 10. iO),lJJ( 3,5,5), CON{iOO),HPi{ 10,10 )» 
i, H F' ... ( ., " , „. I' ) 

Ci:if'lMU’i/TUTAL/,-!CaNST,NOBJ ,NTFN 
-{''^f'CN/Ot'RLH'I/FUMVAL 
CUM'HIR/CYC/CYCLT, AMAXO 

CUMMO.X/PRO/PRQD • 

tU’^HCK/riVAR/I YXF.INYE 

CQ^'YllN/'VaV/AfOMR.COMR ,CCMD, CMC, TCM,COST,TIMEC, TOOL, FORCE , NSPO 
Cun.-lOrj/flVAR/KAMAT.iH^COMR.RSTN.NOPT.NCORD.NCMC.NTCM.NCOST.NTIME.NTP, 
iMFP,fMMACH'j,NSPf)DL,NAXU,NCROSS,NPROB,HN,NSTION,NOP 
COMMCm/irifJ/INUtX,INDEXF,INDEXM 
COYMON/UUT/Y, lOY, iaX,NJ,Gl,G2,FC,HPl,HP2,T,TC 
C0MMC!^/'4A:|‘;/R ,NiN 
CaMMO:j/IN..ULT/CON,SATSFY 
DIMENSION TLIFC.i, > 

DIMENSION TVARUIO ) .FVARUIO) 

C0MM0N7LIFE/TLIF,TVARI ,FVARI,NTA 
COMMQN/TUR/FCC(i . ) ,Gil (1 ),G22{lfi) 

DIMENSION PROBHPIS) 
iSJTA=0 

I F{ INOEX.cQ.O.AMD. IN0EXM.EQ.2) GO TO 25 
I F ( INDEX. :.U.l .ANO. INOEXM.EQ.l) GO TO 31 
I F ( I NDEX. EQ.C. .AND. INDEXM.EQ.I ) GO TO 400 
CALL INTAL(X,Y) 

C INITIALISATION OF X AND Y 

INDEX=C 
GO TO 27 
31 RETURN 
25 D023I = ;.,I0Y 


IDX=IUX+I 

Y{I}=X(IDX} 

23 CONTINUE 

27 IF { INDEXM.EQ.l) GO TO 4L0 
NCONST=NCONST+l 


IOAX1L=0 
NSPNDL = l , 
12 I=NSPNC'L'"" 
NN-NSPC{U: 
D02J==I,NN 
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1 1- { '.u L' ;■ C ( i , J , 9 ) . „„ U . . ) GO TO 2 
I F { Oj 'T ( I , .j , .) (jg TO 3 
1 F ( CliO C ( I , J , o ) . CT, x . ) GO TO 4 

NT,^=-mT 3 + „ 

IF(i.9T.., ) CCFiMl, J,7)=2. 

I F { I . 0,':.: . ',.jn.Caf^D{ I , J , ) .EQ. l. )C0MD ( I » J,6) = 2 . 

T ^iA = l + .i.-I :'j ?, X I L 
PHIH.. L ) 

PUlP ; = Fr,j . 

CMP=T';'M I J J J «« (-5) ) 

r.x=. 

M J ( : , ! , J ) = I NX 

lUX=IuX 

*'lLUi\= iCUN + 1 . 

MAXiKLP F , D CONSTRAINT 
CLiNIi.'CCN) =ALGG(CMC(R!f, 8) )-X(lNX) 

NCn.. = ,\C:JN + 

CiJOCMCC ■■i)=-ALC'G{ C-MC(MN ,7) )+X{ INX ) 

S 0 70 5 

irjx=icx 

NJ { „ ,I , J ) = Ir-iX 

\jTA=NTA+^ 

IF(CUf'C(I,J,7).GT.t.) G0T0 6 
INY=I+1 
I NY = 2 

NJ(2tIf J)=INY 
IOY=INY 
NCrJN=KCUN+ L 

MAXIFU<« SPCcO CONSTRAINT 
CDN(MCCN) =ALCG(CHC(HN,6) >-Y(INY) 

NCDN = lTCON + - 

MIfJIMUK SPFF.D CONSTRAINY 

CL)N(NC0N)=-AL0G( CMC{ RN,5) )+Y( INY) 

GO TO 8 
INY=IOY 
NJ(2tI ,J)=INY 
GO TO 8 
I NX=1 
NTA=NTA+1 
IOAXIL=IOAXIL+1 
MCON=NCOM+.. 

MAXIMUM FSFD CONSTRAINT 
CON(NCCN) =ALOG(CMC(MN, 20 )-X{l) 


NCGN=NCON+i 

MINIMUM FcED CONSTRAINT 
C0N{NCCN)=-AL0G{CMC{MN,i9) }+X(l) 
IFtCQMCd ,J»7}.GT.l.) GO TO 8 
iNY=n-i . 

INY=2 ' ' 

NJ{2t I,J)=INY 
lOY^INY 

NC0N*=NC0N+1 , ..'V'"-:.;: 


WMSk 


CAKE 

CARE 


CARE 


CARE 


CA-e 
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c $!> ''U Ci'iNSTRAI M 

C i..: LOG { CMC { VNt, 6) )-Y { INY ) 

‘vl. ..i JC'/'J + . 

c ;. "jY S'-j.-'-.c Ci:>istrai ky 

Gl./.(\iCC’!)=-ALi:G{CMC(M:'!,5 ) )+Y(INY} 

8 MCG . = <\CG.-]+ , 

G =hyil[l >J,3)Y*Z+{TC0L(I,J,6)*YnNY) ) #=i=2 + ( TGOL { I , J , 2 ) *X ( INX ) )=«'=«' 2+ (TOO L ( I , 
; ( Tl GL ( I , J ,*'■ )«CL1'' J{ I , ) )**2 

G =(TCCL( i, J, n yaLGGC.C. ) )** 2-TCCL( I, J,8)*«2+G1 
C TOOL LIFi; Crj,^iGTRAlNT 

C.iraiVwCiJ) =GL0G( 0. )=i'TQCL(I, J,T> + (1.+TGQL( ) *X ( INX ) + ( 1 .+TCOL ( I 

J , G) ) *Y( riY) +TCQL ( I , J ,3 )*COMD{ I , J,4)-Y( 1 )+PHIPi*(SQRT{Gl) ) 

GLiiK'iCCG) =CON(r]CON)-ALCG(COHD( I, J,3) ) 

TLlFt: = CJi-J(:jCOri)-PHIPl>!=SQRT(Gi)+Y{l) 

TLPFiRTA) =GXP{TLIF2) 

TVkKI {NTA)=TLIF{NTA)*SGRTCG1) 

G;. =FiJRC;. ( i , J, G) **2+(F0PC.{I , J,6)^Y( INY) )=!‘*2+( FORCE! I J Jt21=«'X( INX) )**2 
+ C FORCE ( I , J,4 )*C0MD{ I,,J,4n**2 
G... = G:'-FnRCL{ I,J,c)=!‘»f'2+ (FQRC£(ItJ,8)*AL0G(10. ) )*«2 
MC(.1M= JCOO+'. 

C MAXIYLP FURCc CONSTRAUT 

CO.GCfXCri) =ALGG(TCM( I , J ,8 ) )-ALOG(iO. )X=FORCE{ I , J,7 )-FORCE( I , J,5 )*Y{ IN 
IMYJ-FORCE ( I » J,- )*X( INX )-FORC£{ I , J,3) =i=CQMD( I , J ,4 ) +PHI P2* ( SORT ( G2 ) ) 
FC=(iO.>i'*FORC:HI,J,7) )*(cXP(Y(INY))=«=*F0RCE(IfJ,5) )* { 6XP ( X ( I NX ) ) «* 

IFORCt ( I, J, I) )*(=XP{COMC{ Ii J,4) )=i'=!<FGRCE( I»Jt3) ) 

FVARI (NTA}-FC»!'SaRT(62) 

G22(NTA )=G2 
Gi.'i(NTA) = G: 

FC2(NTA)=FC 

HP„,{ 1 ,J) =CHP’i‘FC*EXP(Y { INY) )*(COMD(I ,J,19)+C0MD( i,J,20) )/2. 

2 CONTI NIE 

\JCON=NCrjN + . 

C SURFACE ACCURACY CONSTRAINT 

+ F (C0I«C ( + I J ,13 ) .EQ.O. ) CON(NCON)=-X( i )+ALOG{ 2. *COMD ( +» J » 16 ) 
l/( SIN{2.YTCM( I, J,5 ) ) ) ) 

IFCCUFUt I, J,i3).5Q.l. 1 C0N(NC0N) = 2.*X(1)+AL0G{C0MD( I,J,16)*8. 

1*TCH( I , J,M) ) 

NSPNDL=NSPNDL+1 
I F(NSPNDL.GT.NSTION) GC TO 11 
GO TO 12 
11 NC0N=KC0N + 1. 

PROBHPd )=0. 

ShP=0. 

D040iI=*i,NSTI0N 


40 

401 


MN=NSPC(I) 

D040J=l,NN 

I FICOMCd , J»9).EQ.O.} GO TO 40 
PROBHP (1) aPROBHPd )+ALCGC10. )«HPX( I, J ) 
PROBHP (1) * ( PRCBHP( I J*FCRCE( I» J»8 J 
CONTINUE ^ G.. ■ 

CONTIN’UE G; G'Y 

SHP = SHP+PR0BHP'{'3.J'*PHIP2,: 

PR0BHP{2I«9.'f: 


.V * . 




ii (. fKif. inji-v 







DiJ''.- I = . jMSTIC i 
N. = .viPC(n 
Oil--.. J=*,n,! 

I F (CQy,:( i , j , -j) , ) Q.- 70 4^^ 

PRfjiHF ( j; ) =PRuv.HP { i ) + Y( K) *HP* <1 » J ) 

PRl>.,H?( ) ={PR'.,RHP{ i>*FCRC£( I,J,6 ) )=!=*2 

41 COFlTiMI-r 

4i: cO'.timj:: 

SHP = SHP+PM;jFHP {.. )*PHI FI 
PP.UhHP{ 

Dl',!4'...I= ,1-1 ST I ON 
'■‘JN = .\lSPC (I ) 

DU4_ J = 1 t'TI 

iF{ccyij{ I, j,9) .:q.o. ) gc td 42 

K = !‘jj( i,l , J) 

PROBHP( 3) =PRn4HP(3 )+X{ K)4HP1{ I, J ) 

PPfiaHP (3) = (PP.JBHP{ ..)*FCRCE{ I, J,2 ) )«*2 

42 CDRir'it"' 

422 CONTIKU:. 

SHP=SHF + PRUBHP( P2 

PROBHP('':.)= . 

004331=} , NSTION 
.M(^ = OSPO{I) 

D043J=ltN0 

IF {Ca^D( !, J,9) .EQ.O. ) GO TO 43 
PR0BHP(4 ) = PR0Bi iP{4)+ COPD( I i J » 4 ) *HPI C I , J ) 
PR08}-IP(4) = { PROBHP(h)*FCRCE( I, J,4 ) )**2 

43 CONTINL't 
433 CONTINUE 

SHP=SHP+PRDBHP(4)4PHIP2 
SHP = SyRT( SI1P)*SCRT(PHIF2} 

CQNJNCCN) =CMC (HN,X1 )*CPC{MN,15 r + SHP 
D09I=*,NSTI0N 
N = NSPO(i) 

DOl.; J=1,N 

IF (CCMCII , J,9}.EQ. ... ) GO TO 10 
C HORCE POWER CGNSTRAINT 

C ON { NCON ) =CON { NCOH ) -H PI ( I , J ) 

10 CONTINUE 
9 CONTINUE 

D026I=1,IOY 
I DX=IOX+I 
X( IOX)=Y( 1 ) 

26 CONTINUE 

FUMVAL=FUNXON(X) 

DD45I=1,NSTI0N 
NN = NSPC(n 
D046J = 1,NN 

IFICO^DI I , J,93 .EQ.O. I .GO TO 46 
4 NjCON=NCON + l ■■ 

CON { NCON )=CYCLT-T( If O' K*01 < 

46 CONTINUE .... ^ ■}■' 



ooo ooooo 
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45 COnilML 

c '‘-i'-'* C'ni> 0 \,,jTs CnNS7R-INT 

ZUiU (CC; i) = ‘.LC ”i{ COMRt i ) )-Y{ I. ) 

;‘s|=I‘J:=;+Imy 

RIzTUR^'i 

t ' ‘I L) 


suBRCtir - i : MLMif -* A ( x ) 

C; ij*: i*c ::Js sjj? ^ ^ '^s 3js ^ 3 *c ^ jjc s^c 3 ^ ^}t :ijc :!:^ ;^ 

C * THIS SUlUCUTI.Jt IS THE BACKBONE OF THE PROGRAMME . IT CALLS VARIOUS * 

C * SLBROUTIfii.S ^: Tf ) FUN .SUBPROGRAMME D IRSCTLY / IND IRECTL Y * 

Q ❖’}= ^ 4: >i* ^ T* ^ >i! ^ ;jc ^ ^ :5c 4 4; ^ jjc::^! :J: #sjt4::!jcf»! 4:* lit# si 

CfjMMCN/DI V/Cl 
: COMMOi /- 1AHE / N,NiN 

CGMMUN /' 1Iia/R 
coMMO ' M/Pov : l/z 
cummon / serche/funval 

COMMLiH/IN::QLT/CnN, SATSFY 
DIMcNSIUN CON{„ > , X( 5C ) ,Xj. ( i.jv ) 

£ :{c^t :<c>ic: 4 i:{t :}c 4 :ii* 4 t:ic 4 ;:i!t cjecjccjt:}!^: cic ci: :ii :}:it cjt 4!:^:it:|c:!!^e:i( ^;*:ii:ie:{c:it:{:4::!t:i!:J;4::5:*:i:* ^::^:{c:{:*:{t :it#:ic:{c:}!:ie* #:i£##*:Jc :4!:J::^; ^ 

C * THF : initial VALUE OF ^ R * IN PENALTY FUNCTION PROGRAMME IS TAKEN 10 >«= 

Q :i(4e:^:4c^3icci'^C)<4ccit:i::{t4c:ic4^cicc^c!‘^^^4>ic:ii:i::^:;c;:!(:it:it:i::^:5i:St:it:ic:^:i::^:ii :i::ic^:ic::ic:i::i::i::ic:ic:^ :{c cjc 4< 4 ^ >ic sjc cic 4 : sj: $ lit ;ic cjt 4c:i: cjc : 5 c :}! ;ie 

I=t 

CALL PINT { X , Z ) 

IF { I . NE.;j ) GUTO 53 
PRINT lOOiFUNVAL 
CALL OLTPUT ( X ) 

CALL CHECK ( X ) 

IFISATSFY -. EQ . .) GOTO 53 
PRINT 59 

59 FURMATf //^ X,#INITIAL SOLUTION DOES NOT LIE IN CONSTRAINT SET *) 

RETURN 

53 CALL PCWELLIX ) 

I =1 + i 

IFd . EC . l ) GOTO 89 
DIFF = ABS { Z-Zl ) 

:Jt:ic***#i^#*$:5;:it:i(:5«:itt*«****^* *************************** *********************** 

* THE TERMINATION CRITERIA 

* TERMINATE WHCN DIFF. BETWEEN PRIMARY AND SECONDARY FUNCTION BECOMES * 

* LESS THAN 4 . U**-6 TIMES THE PRIMARY FUNCTIONAL VALUE * 

*«#:5c^!*^:^*:5;#*:i!* *********=»************************* ************************** 

IF { OlFF . Lt . ABS ( Z * l . E - 03 ) ) GO TO 9 © 


80 Zi=Z 

♦ 4; ***:{«:)! iic:ic:it!ic«*:5!!ic*:ijc#:«::5: #«*««*4:**$*<c*«:^*«*:^**«*«****«*«*****=4«* *;»=»* 

* PENALTY CONSTANT *R* REDUCED TO 1/4 OF ITS ORIGINAL VALUE * 

♦**:{t«:{c:i(:)(*:ie:i::{c:5c##«4c!je«#*:!c*:}4************ ************************ *************** 
' R=R/Cl ’ c.' ,, . ^ , 




' is‘4 “ B , 

.-'S ;■■•■ ''j 

.Tfi ' B '-- -1 



, , 'i:' "'n' J 

■* . 



'''fO'i VV!"'' ‘ 





o o o 



FORMAT (///, .. 
.,R «) 

R = t : . 


2 ', 

\,*IMITI4L FUNCTION VALUE*, F15. 6, 5;:X,*SOLUTION V£CTO 


CALL PCU;:LL(X) 

PRI‘.T /5 

FuRR.^T ( //////// I , X f THc uPTIMUf' SOLUTION HAS REACHED 
PHIRT 2 ,(CriU( I) ,l = ;i,NIN) 

FORMAT C.X, OF 0.4) 

PRINT 5 tFUNVAL 

FORRATI .'(j^orTIRAL VALU:- QF OBJECTIVE FUNCTION=*, FI'. .4 ) 
CALL CLTPUK.O 
Rr. TURN 


S'jBRQUTiN., GR^D(X,&) 

£ ' sj: ^ ^ ^ jjt: 5 ^ 5 ^, j§c s|s 3 |s sj: : 5 c : 5 c 3 {e ) 5 < : 5 : j^c 3 {c' 3 ^ : 5 c . 3 ^ ; 5 « s^c " 55 ^ : 5 c 

C *TillS SUBROUTIN'C EVALUATES THE GRADIENT BY FORWARD DIFF. FORMULAE * 

Q ❖ 4 ********»r**=;' !{!^s^^ 3 ic 3 js sjt:^ jjc ;{t :{e ;!t ^ 3 }: 4 :: 4 ! :{t * 

C0MM0N/MA1R./N,NIN 

COMNCN/IN;.QLT/CON, SATSFY 

COMMON/PUVEL/Z 

DIF.CNSLTN X ( 5 0 , XI ( 50 ) ,G ( 50 ) 

DIMtNSIQN CON{^ .. ) 

FX = Z • 

DO if;-,, 1=,.. ,N 
DO iO J=1 ,N 
10 Xi ( J)=X(J) 

DX=X( I )*J.. 

IFIUX.EQ.O.) DX=0.01 
5 XKi )=X{I ) + DX 
CALLCHECKIXl ) 

IFtSATSFY.Cy. . ) GOTOZQ 
IFIDX.LT.O.) GOTOll 
DX=-UX 
GOTO 5 

il DX=-DX*R.l 
30T05 

20 CALL FLNT(Xl,FXi) 

100 G(n = (FXl-FX)/DX 
RETURN 
END 


FUNCTION FUNXONIX) 

«* 4 :* 4 <«*:»>)!*«***«**«:}:** 4 *=t* ***********************'*'******♦******************* 

* THIS FUNCTION SUB PROGRAMM EVALUATES THE OBJECTIVE FUNCTION * 

:}:;}«**« 4 :* :«:*!«!#********* = 4 **************************=^*******^**************** 

COMMON/TOTAL/NCONSTtNOej,NTFN 
C0MM0N/CYC/CYCLT,AMAX0 ,., ; V 

CQMMON/MANE/N*NlNM‘"''''^'‘T"-T-'';" > '-N' ; , , . ; 

C OMMON/MATV/AMATR»CO«ftfC:dMpjCWCtTCMt COST* TIMECf TOOL t FORCE, NSPO 
COMMON/NVAR/NAHATlKW^ptNtNOPT.,^ ^ ' ; 



o o o 


J >MFP,,:vicj Ij, ^S|)^^rjL,NAXI L,NCR0SS,NPR0B,f#4,l\iSTI0N,iS10P 
Ci.i*'-'': /I :'4D/ii'.;D:‘. I NDEXF, INDEXM 

C'. ''''-'/'^*J!/X,i.i}YflOXjf'siJ,Gl»G2»FCrHPl»HP2fTjTC 
.;/!rr-,;-:LT/CCN,SATS FY 

D IIJM 4SPC{ . u) , AiMATRIIO ) ,COMR(lO) »COMD { 8 , 6 , 25 ) » C MC ( 5 , 2‘.. ) . 

:TC",(;-,6, ,CQST(] ) r T IM :C ( 8 , 6, 1< J , TCCL { 8 , 6, i 5 ) , FORC £ ( 3 , 6 1 1 5 3 t 

),y(5 ),fp*(*>:,i ),HP2{lv',ia),x{5v. )tT(i^,iw;),TC(ic ,ic) 

t * R 2 5 ) 

TY--.T (50 ) ,TMATC{50) 

I F { I !\OiX.-i. -u.., ) RETURN 
■rLiHJ = \C-’.J + „ 
is:TFG = .-;TFO + ,., 

CYCLT.=.Z5 
Kil = 

C05T( X)= , 

L* 0 ji. i ~ ^ » i '1 S ! I Q !’ii 

*^iM='lSPG^I ) 

D u2 J = i t 'JN 

IF(C0PC{I,J,9).EQ. GO TO 2 
INX=riJ („, I , J 3 
INY=NJ(2, 1 , J) 

T ( 1 , J 3=C0i'1D(I , J,3)/(EXF(X(INX) )*EXP(Y{ INY3 ) ) 

COSK ;i0 3 = CrjST{l )+T(I, J)^COST(l) 

TC ( I , J )=AbS(T(I , J)-CYCLTi) 

KN=KN+I 

TMAKKN 3 = T( I , J) 

TMATC{ KM)=TC( I, J) 

2 CONTINLE . 

I CONTI NL2 

CALL MAXH{TMAT,CYCLT,J ,KN) 

GO T0(6no,700t800) ,INDEXF , 

600 FUNXON=CYCLT 
RcTURN 

7 It: CALL CBJFN(TCOST) 

FUNXON=TCUST 

RcTURN 

6W CALL CBJFN(TCOST) 

FUMXON=TCOST 

RETURN 

END 


SUBROUTINE POWtLLIX) 

=i' THIS subroutine PERFORMS 0-F-P OPERATIONS * 

sic sis*** *****=*>*****=> >!'*************’^*^******=*'****=^* 


COMMON/MI NI/R 
COMMON/PUVEL/Z 
C0MM0K/SERCH2/FUNVAL 

COWMON/S6RCH3/F ■ ‘ /, , , , 

COMMON/MANE/N,NIN ^ 

OIM£NSION'"x'{5^!?*6l5oI!m20f2OI*SaO»»SIGMA(2OI»YY(2OifHl(2Cf,20), , 





O O o ' o o 


c '• ) 5 ' j 50 

0. •■' ..'I : ) 


c * Ti-: ..^:L':. ..'.ivis "hi accu acv cf step NCRi^i/'n 

15 . Pil" ' CO--. ■■ .*-M- O 

i.i'V .. I-., , , 

'j . .. J I ■ ■ 

■i I i j J ) = • 
s: HO , I ! = . . 

C/-Ll. .,;P'0(' jG! 

1 1 \ = 

9i Gul OO 

U Ui t 1 ^ ♦ 

Li — , # ' , ' ’ * 


0 j G. I =■ o 

S ( 1 ) = . 

DC ) . J= 

91 3(i )=5{i )-!'( I, 

[0-i;+'i { I )>!'3 (!) 


IF {U.Gi: . . 5 O: 'Ti., 7. 

9_ UijT9S = S { I )=5«S { 1 )+L'lTSS 
DC 99C 1= 

S( I ) =S ( I ) /SC'^T ( CCTSS) 

99^ CQiOIMv, 

CALL SEPRCHCOo, UFI ) 

COI'JV .PCcnC. CPITCPlii ST5F SIZE 

IF(ABS(^LP^ .LT.(’.j.E-J2) 9LTURN 
C VAL'Ji CF ALPH,fl = . IN LINEAR FIN. FtANS NC PRfDGRESS IN LINcA 
I F( rfN.;;U, ) GCTC 91 ■ 

IF (Z.LT .Z ) GLTC 9,. ' 

C THIS FCLLOHv F ( ^(U + AL PHA>!= S ) = FO 

z = z,. 

RcTLRN 

92( DO 9 2 I=i,K' . 

SIGFA { I ) = &LPHS{ I ) 

CT = AMA5<1 { , APS { ALPI ) ) 

DO 75 L=i,N 
C 

C CUNVERG£NCF FCR lOTAL. F0V5 
C 

IF ( ABS {f:T=^S(L) ) .LE.FPSLCNd) ) GOTO 76 
GOTO es 
76 CONTI NLc 

PRINT SCJ'.itRj ITN,Zf F,ALPI» { SC'J) » J=ltNI 
CALL CLTPUTCX) 

RETURN .. . 

89 DCTSIG = OOTSIG+SIGt»!A( n =*’•^2 - ■ - - -V. 

92 XdlaXdl+SIGFACI) ' L''''-"'.';,: 


R PIN. 



LL 


r ’ = •: -MG { 3 !*'i ( j ) 

Y Y ! . ) = G I . 5 - 1 . ) 

■. 3 ( I) * . ( . ) 4 

9': '3 c ; ) = . : (1 ) 

C 

c ’>i.. .:UG Gl.iT- "4 

c 

3 r {i-::/ n 
A = , 


Nl IS TO SKIP THE UPDATING OF H-'YaTR-I^ IF Pi 
'^-0 A ) GO Tc ee 


fi A= . ■ " , ^ 

D-: '^5 i=.„ , ’ 

a: =..., +sk;> ^(i )=:=\y( ii 

Dl'i 95 0=.?'. 

95 3 . -..,.-, + 'irY(!)*hl. ,J)^'YYf J) 

t .= ,/ a. 

no 9;i 

D ! ) 9 J ) 0 

II ,(1 ,J)=H{I,J)+3lGM,' (1 )-'!=AIG^'A(J)^vAi 
DU 9c K=., ,!3 
DU 9e N = ., jO 

98 H., (I » J 3=H, { I , J ) +A?*H( i ,K 5 =^YYI K)^YY{ F J) 

DC 99 I =1,0 ■■ 

Du 99 J=' , N 
99 H{I,J5=ri3. (I ,J) 

3 6 I 'i ! ' - i 1 i + ; ■ ' ^ 

1 . ! ■( I 1 I' » ' i‘ I / , 1 , /* L f' i [ A ( i } s , F ) 

• t r,». ^T'l ] " ! '• ) 

TF(ITK.la.-YN) GO TC lU6 

C THUS 2C IS THf. ^AXY, NC. OF ITERATiONS FOR A VALUE OF *R=!= 

ICC Z.=Z 

GOTO 9C 
H3 6 PRINT 1 7 

IU7 FORPAT (///J CX ,*^AXIFUW NC OF ITERATIONS CXCcLDLC HuNCt QLIT*) 
5CC0 FORMAT { IX , = * , r-^2 . 6 , 1 5 , l^Fi 1 .4, U F7.3 ) 

i<35 RETURN 
END 




SUBRQLTINl IMT(X) 

C 

C THIS SUBROUTINE TO INITIALIZE THE VARIABLES IN CASE 0 F SINGLt- TCUL CASfc 
C 

dimension X(5s5 
X{i)=2C0. 

X{2)=.02 

X{3)=J30. ^ ...V,, ,, . , u ■ - . 

X{4) = .2 ' ■ 'u ;V'Y, 

' , . hhturn u",/' ■■ ,v ' 












' r' , 

c/' 







--v ;! .'.yr ■ c,>, , 


T'U 9: • 

‘ v’/". 






' ■i'' 





■■t 

; ■: ' ■ .f/'i 


: 



o o o o 




i. T * i 

ri..i(,.,rx) 

c 

V ;**» "1 

% 'i */ '1' 

.•!■ ;I< ■■? -ii a a* 4 4 4*4' 4 4“ 4 4=4 4= 4 4 4 *4# 4*45 *4: 4= 4= 4 4 4 4=4 4=4= 4 4 4* 4i 4: 4;4 4: 4^ 4 4= 4 4 4 4 4 4 

c 

j t" 

S > 

^ the sccgmcary functicn 

c 

<lji '■»* «•*' 

'.Ji ':J; [c 

;■ 4 .4 4 4 4 > 4 ;;< 44 ,-, x 4 44 44 44 44444444:4444 444 4445lt4:44 4444=44 444 4444<4!44 4 


L .’v. 

>./’■: 

/i^ ; ■ 


) V */ ' 

■,/■■■' ' 

/C , iV 1 


r ! ■ ' 

N/: i 

. '.'Ll /-..GY, SITSFY 


C .. L 

0/ ... 

a.h ;;/FMNV.iL 



I""* 

I'-i . 5:0 j :':a!. . 



oar {. ) 


F 

,‘'=f-L'-V .L 


Tj 

: ! = ,. ,Ji!4 

iC 

F 

J = F.\ + R/C.n ( I ) 


R 

Tl'Lv 



, 1 r> 

i L,.i 


S!J‘iKCLTl Ch .CK (X ) 

^:{£ 4 jj! ^ 4 ' + -l! >!! 5 !< 4 j!' J s!' 5 ;t 4 =4 ^ ^ ^ 4 J{: * 4 4 >!< 4 it ^ 4 5 !=# V ^ ^ ^ 4 ^ ; 

* fl.iis St.’RClTIXH iS 10 DHcCK ThE FE^SIfilLiTY Rt:GlQ;>i 

4 < 4 £ij!)l[ 4 ' 4 ' 4 i 4 iI‘ 4 ‘ii*TT'ir-TT''rii‘Tiii' 4 'iT=l‘>!i 4 'vi?Tii 4 = 4 < 4 s 4 'ii=i 4 i 4 i!‘i 4 i 4 iif 4 'i 4 i 4 = 4 ': 4 * 4 =i!'i!ii!> 4 ' 4 ‘ 4 ‘i!' 4 ii 4 ii{'ifi^i!'iri 4 =T'^ii‘ 4 ii 4 ‘ii' 4 iY 4 =: 

cnf^r'tN/ 4A .. /R IR 
CC^'YCN/ll'iEiLT/Cofv.S.jTSFY 
DI H-.iYS lO'j Cf.:K(.. JO) 

Oiy, ^ 5 Ij;,j ;/(!, ) 

SAT!5FY= . -' 

call CC.1STR(X ) 

DC. .C !=: ,NI L : " V . 

IF (COIN (I). LI. . ) GOTO 15 
iv COOTlNL'l 
RETURN 

i5 SATSFY=.. 

R 'TURN 

EJJD 


SUBRCLTINr S'iARCHIXrS, ENIC) 

4 4" 4^ 4= ij! ij^si' X * iJ !> * * « =«• 4 4=4* 4 ♦ 4^ i* 4=:4i * i?'*^* =«i4< 4= 4; 

4« THIS SUBRCLTINE PcRFGRMS THE GCLDEN SEARCH OPtRATIONS 

!{tjle 4 : 4 c 4 ! 4 i 4 4 5 l>,)t!ji>!<>!t 4 i}!i 4 i 4 tJ 4 ! 5 !£ 45 f 3 !‘!{'){i 4 > 444 ! 4 i#i!! 4 !Ji‘ 44 ii> 44 ' 44 < 4 i 4 * 4 <iii 4 ! 4 < 4 ! 4 si!!* 4 = 4 i 4 i 4 * 4 i 4 '*i 4 !i 4 :ii««i^« 4 : 4 '* 4 i 

C0HHGN/Sr.RCh2/FLNVAL 

CQYYCN/PUVl-L/Z 

C0PiVCK/S£RCH3/F 

COMfYQU/MAMF/N,NU>i 

CnyyON/INHQLT/CQNtSATSFY 

CGKFCN/STtP/ETA 

DIMENSION CnNQCD) 

DIMENSION X (50.0 ,'S{5C ) ■ . 

DIMENSION YY( 5a) 

CONST '1=':>.5R 1966 ^ '■ ' ' , ■ 



I 




r. / 

, , U 

■ ■ I = 5 . 


1 i 

■'Y ; ; ! = (;.)+.■ ,V. 

C-.L'. '> CY(VVJ 

S{ I ) 


3X- -Y . k ) 

j = : 

: ■ ’■ X / . 

T' • 

yj C 7 C 

9'.' 

ufJ.":. ) X.' cc 

..1 . .=t T 

*f - .. 

LJ ; i..,: ... „ 

'7-T 

9X 

Y 1 = . T - 


- ■! 

i ». - 

i.( I- . . 4 i — f \ 


li* 

YV( 1 3 = '{ i) + - ' 

CM.i C:l .(:< { Y Y) 

- S ( I ) 


iFc-^ATcrv. .) 

,, ft . 0 =;•■ : M -i 

GL TC 1 

aC T 

lit 

F !^h'/ = E 


1C 

uisr='.x.4/.- YUs 
OY=i,aST 


n 

2., = ;; ^‘'I^+c,:).•JST 



FX = f,f^I^+C.J•^!5T.lC 



Call f-Lii . ( Xi S j . ,F ^ ) 

FIJ', = PLNVAL 

CALL Fix: ( <,S| ,Fii 

FU2=FU^VAL 

Ic I F{AB.3 {'..MAX-'Iwir.), Lc.O .C'..'1‘ABS(CIST) ) GOTO 3C 
Dc =const:.:-h;:'-, : a 

IF{F,i-F2) 

14 5MIX = 21 . ; X / X ' : 

X]=::J ■ ■ 

F3 =F2' 

FUI = F!J2 

E2 = Ef^IN+C01ST:*Cc 

Call fln: { x,StC.:,F2,) 

FU2=FL^VAL 
GCiTC .,2 
13 Ef^iAX = c2 
F 2 = F ik 
E ? = •“: I 
FU2=FU1 

Fi.=cWI^+ca^iST: *c£ 

Call FLMl (X, S tEXt F1 ) ■ ... 

FLl=FLhVAL 
GUTC :.-2 ■ 

25 1 F( ICC.EQ.l ) GC TG 13 . , 

£HIH = 21 ; . , ‘ 

EI'A'A-E2 . -X'' ■' X' ' 'X-. 





3C irCF. 

^ I [ . = 

Z = F 
r " F . j . 

1 1'. F = ", c n + ' 

1 ^ ( If..,'. r . ) 

. F J N = ’ . 

-f, ; 

K ' ■ I ! = . r 
F = rL : 

ir-nci. .FT. ) 
r Fi r = : . 

[- V C. ; = ^ ; j , 

0 c I F ' .. : 

33 k.TLFiX 


StJl'i3CLlIN: FL^ (>■ , S , c. , F ) 

ct i«ycN/iN,- jlt/clkt sarsFY 

IJIM'-N!: IIJM X( St ) ,S( 30 ) , YY{ 5 G) 
dim,. KSiOls CCKi.^,, ) 

DO 

io YY(i )=X{i) + :;*s(i) 

CALL CCIJSTR{YY) 

CALL riNT(YY,F) 

RLTURK 

EFiC 

SUORCLTII'L: MAXM{t, AYAX ,J,K) 
DlRl:KSl!Jfj ?{*; ) 

AMAX = -10. >l'*4 
I = I ,N 

IF { e ( I ) .Lr.ARAX) GO TO 1 . 
AyAX = E I!) 

J = I 

I.C CCiNTINLL 
RETURK 
END 




